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INTRODUCTION 


The computer program to be described is entitled "Advanced Investi- 
gation of Thermal Protection Systems for Atmospheric Entry-Modification 
Number 2" (CHAP III) . CHAP III is programed in FORTRAN language and is used 
to compute the transient response of a thermal protection material to a 
prescribed heat input at the surface. The original version of the program, 
CHAP I, was developed (ref. 1) to deal with ablation materials but has been 
applied to a wide class of thermal protection problems. CHAP III has been 
developed primarily to add the capability to analyze chemical kinetics in 
detail. 

The purpose of this note is to provide a user's manual for the revised 
computer program. Use of CHAP III should not be attempted without 
familiarity with the mathematical development of CHAP I. Lengthy 
mathematical developments are avoided since they are included in the 
references. However, this manual does provide a brief description of the 
mathematical model around which CHAP III is built. The emphases herein is 
on those points that represent changes in the mathematical model that was 
the basis for CHAP I. 

This manual contains a set of tables which provide: (1) a definition 

of physical symbols and a list of FORTRAN symbols, (2) a list of input and 
output FORTRAN variables in the order they appear in the I/O lists of CHAP 
III and, (3) an alphabetical list of FORTRAN symbols and their location in 
(1) and (2) . Cross reference features in these tables should permit 
straight-forward use of CHAP III. 


SYMBOLS 


The computer program is normally used with English Units 
shown in parentheses below. A consistent set of SI units can be 
used if PREHE = 2 and YDEL is set equal to the Stefan-Boltzmann 
constant in the units chosen. 


frequency factors in Arrhenius constants (units vary with 
reaction order) 


chemical specie identification symbols for ret^tants 


solid. 


chemical specie identification symbol for products 

3 3 

concentration of a component, gm-moles/cm (gm-moles/cm ) 
mass fraction converted to carbon by pyrolysis 


E 

h 

k. 

L 

M 

m 

m 

n 

w 

N 

R 


w 


q c,net 

S 


activation energy, J/kg (Btu/lbm) 
enthalpy, J/kg (Btu/lbm) 

leactlon rate constant tor pyrolyaia reactions or for oxidation 
of char (units vary with order) 

total thickness, m (ft) 

molecular weight 

number of chemical reactions 

2 2 

mass flow rate, kg/m s (lbm/ft s) 
reaction order 

3 3 

rate of deposition kg/m s (lbm/ft s) 
number of reacting specie 
universal gas constant, J/kg-K (Btu/lbm • 
pressure at the wall, atmospheres 

rate of production of chemical specie by chemical reaction, 

3 3 

kg/m s (gm/cm s) 

2 2 

net aerodynamic heating rate to the surface, W/m (Btu/ft s) 
power on the temperature in the rate equation 


t 

T 

V 


a 


ij 


0 

<5 


K 


time, s 

temperature, K (°k) 

velocity of surface due to surface recession, m/s (ft/s) 

stoichiometric coefficients of reactants 

stoichiometric coefficients of products 

order of reaction for reactants 

order of reaction for products 

2 2 

char permeability, m (ft ) 

2 2 . 

viscosity coefficient, N-s/m (lbf s/ft ) 


2 


§ dimensionless coordinate normal to surface 

X weight of char removed per unit weight of oxygen diffusing to 

the surface 

3 3 

p density, kg/m (lbm/ft ) 

* b l parameter equal to _ m 2 a/[| (£t 2 s/lbf) 

ZRTy 

*p parameter equal to m 2 g/t| (£t 2 g/lb£) 

ZRT|i 


Subscripts : 
c char 

e external to the boundary layer 

f forward 

i integers 

o initial values 

p pyrolysis or constant pressure 

r reverse 

s surface 

w wall 

MATHEMATICAL ANALYSIS 


The Model 

The physical model of an ablation material as used herein is 
shown schematically in figure 1. More complex models have been described 
(see, fcr example, ref. 2). However, the physical model used herein can 
approximate a wide range of ablation problems. Moreover, this model is 
described adequately by a mathematical model whose computer solution Is 
feasible. The general form of the energy equation for CHAP III is shown 
below: 


(Storage) + (Pyrolysis solid reactions) + (Pyrolysis gas reactions) + 
(Convection) + (Conduction) + (Work) + (Kinetic energy) + (Viscous (1) 
dissipation) * 0 

The mathematical model differs in some respects fro* CHAP I. The 
areas of the analysis where major changes have been made in the terms of 
the energy equation are outlined in the following sections. In addition, 
a set of density equations has been added to account for a pyrolysis 
reaction zone rather than the reaction plane used in CHAP I.- 


The Pyrolysis Gas Kinetic Reactions 

CHAP III can compute the solution to the generalized chemical 
kinetics problem. The kinetics analysis is necessary when neither 
chemical equilibrium nor chemically frozen pyrolysis gas flow provides an 
adequate approximation. In addition, simplified analyses are included to 
provide the bounding conditions : gases that are chemically frozen or 

gases that are in chemical equilibrium. 

The general kinetic solution .- The kinetic solution must be 
approached from a general standpoint to have broad applicability. The 
chemical reactions involved in a generalized reaction can be written in 
the following form: 


N 

Z 

i=l 




N 

Z 

i=l 


6 ij B j 


( 2 ) 


Each value of i in equation (2) corresponds to one chemical specie and 
j takes on values from 1 to the number of chemical reactions involved. 
Once the equation set. has been established by equation (2), the rate of 
production of the i specie can be written as: 
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P 
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- Z 
J“1 



V (k f. 





(3a) 


where the values k^ and k r are of the form 


k 


AT 


-S -E/RT 
e 


(3b) 


for each chemical reaction. 

The generality of the above equations makes them very useful. The 
price paid for this generality is increased computing time. The chain 
multiplications involved in the it terms of equation (3) are the main 
contributors to large computing times. 
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The above equations are included as a sub-routine that can be 
inserted into CHAP III to model a generalized set of chemical reactions. 
However, the generalized approach is not recommended except under unusual 
conditions. 

Explicit kinetic solutions .— A more practical approach to the 
kinetic solution is to write out explicitly the rate of production terms 
defined by equation (3) . This approach constrains the set of chemical 
reactions that can be modeled. This constraint will be minimal if 
judicious choices are made for the set of chemical reactions. For 
example, a system of chemical reactions for ablation materials has been 
programed explicitly and used with excellent results by the authors for a 
number of ablation materials. This explicit equation set is included in 
the program as the standard option for modeling gas kinetics. 

Simple systems should almost always be modeled by writing a new 
explicit expression to describe the unique system being modeled. This 
explicit expression is programed and inserted at statement number 10 in 
the KINET subroutine. The appropriate changes can be made in the comment 
cards in the INOUT subroutine so they describe the set of chemical 
reactions currently being used. An example is shown below to demonstrate 
the simplicity of the explicit approach. A simple set of chemical 
reactions is shown in equation (4) 

^ CH 3 + H (4a) 

^2 ^ 2H (4b) 

The 4 species are: 

1 - CH. 

4 


2 - CH 3 



4 - H 


From equation (3a), an explicit expression for the rate of production of 
specie 1 (CH^,) contains only one term since CH, appears only in 
equation (4a). Using FORTRAN symbols: * 


Y Y 

XMULTP - C 4. C 2 

H 3 1 « C(4)**NPEX(4,1)*C(2)**NPEX(2,1) 
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XMULTR « C 


C (1) **NREX(1 , 1) 


( 6 ) 



p x = KATE (1) “ < p SC (1,1) " RSC(1,1))*(FK(1)*XMULTP - RK(1)*XMULTR) (7) 

The expression for RATE(2) is equal in magnitude to RATE (1) since 

CH 3 appears only in equation (4a) and the stoichiometric coefficients are 

equal. This is programed as: 


? 2 = RATE (2) = -RATE (1) 

The expression for RATE (3) also contains only one set of terms. 
Y 4 

XMULTP = C R 2 = C(4)**NPEX(4,2) 


XMULTR = 



- C(3)**NREX(3,2) 


( 10 ) 


P 3 = RATE (3) = (PSC(2,3) - RSC(2,3))*(FK(2)*XMULTP-RK(2)*XMULTR) (11) 


The production rate for specie 4 contains two sets of terms since it 
appears in both equations (4a) and (4b). 

RATE (4) ■ (PSC(2,4) -RSC (2,4))*(FK(2)*XMULTP-RK(2)*XMULTR) 

+(PSC(1,4) - RSC(1, 4))*(FK1)*XMULTP-RK(1)*XMULTR) 

It should be noted that the expression for specie 4 can also be written 
as: 

P 4 - RATE (4) - -RATE(l) - 2*RATE(2) (1 

If equations (5), (6), (7), (8), (9), (10), (11), and (14) are used, the 
computer time is approximately 6 % of that required by the generalized 
equation (3a). Larger systems of chemical reactions can be programed in 
an equally simple fashion. The approach to a kinetic solution in this 
section has followed closely the development in reference 3* 

The boun ding solutions .- It is frequently desirable to establish 
the boundaries within which a solution will lie. With respect to 
chemical reactions, these boundaries are defined by: (1) the solution 

when the pyrolysis gases are chemically frozen and (2) the solution when 


( 12 ) 

(13) 
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the pyrolysis gases are in chemical equilibrium. The frozen solution 
requires no elaboration. The equilibrium solution contained in CHAP III 
is based on the assumption that the pyrolysis gases are in local thermal 
and chemical equilibrium within the porous solid. This assumption makes 
it possible to compute the enthalpy and chemical composition of the gas 
as a function of temperature and pressure. The equilibrium solution for 
a given chemical systen is readily obtained by use of the computer 
program described in reference 4. The output of that program is 
compatible with the input required by CHAP III. 


Chemical reactions between pyrolysis gases and boundary layer gases. - 
CHAP III hat, provisions for incorporating the effects of boundary layer 
combustion into the surface recession rates. The :aodel used to develop 
the combustion analysis was a carbon bearing gas being injected into a 
boundary layer that contains oxygen. When the combustion of pyrolysis 
gases is included, the expression for first order oxidation of the 
surface becomes: 


m 

c 


Here ® car jj 0n includes all carbon injected into the boundary layer, 
either as a result of pyrolysis or as a result of sublimation, except 
carbon in the form CO. Also, the assumption is made that any oxygen in 
the pyrolysis gases is in the form CO. The computer program has been 
made flexible enough so that other material systems, exposed to more 
generalized atmopsheres , can be analyzed. 


KP C - 
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(15) 
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The Density Equation 


The chemically reacting components of the solid are allowed to vary 
as a function of position and time in CHAP III. Pyrolysis reaction- in-depth 
has been incorporated into the transient analysis in a fashion similar to 
that used in the quasi-steady solution described ir reference 5. This 
analysis is not included in reference 1, so it is outlined here. The 
differential equation for the densities of the Irreversible pyrolysis 
reactions are: 



Dt 


dp 

—1 + 
3t 



a§ 

3t 


(16a) 
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where there Is one density equation for each component in the solid. The 
rate of mass loss 3p/3t of each component is jiven by a pseudo-order 
classical rate expression (refs. 6* 7). 


9£ 

at 


n -AE/RT 
p Ae 


(16b) 


Solid carbon is formed both by the irreversible pyrolysis reactions and 
by deposition from the gas phase. The differential equation for the 
density of carbon is: 


DP. 


Dt 


NTYR 

Z 

i=l 



^solid^i 


+ .w + 



3§ 

at 


(16c) 


The rate of deposition of solid, w, is calculated from the gas phase 
kinetics solution. The expression for the derivative 3J[ is: 

at 

If ■ — a - i> (i? > 

dt p L 
s 

where m is the surface mass loss rate and p L is the local density 

S 8 

at the current surface location. The remaining derivative in equation 
(16), j)p , is obtained with standard finite difference techniques. 

a§ 


NUMERICAL ANALYSIS 


The physical model in figure 1 has been approximated in 
one-dimension by a finite difference network. A coarse grid illustrated 
in figure 2 is used with an implicit-finite difference formulation of the 
energy equation to generate a set of simultaneous algebraic equations. 
These algebraic equations are readily transformed to a tri-diagonal form 
whose solution is rapid and accurate. Within the coarse grid of the 
energy equation is a finer grid whose spacings vary with chemical 
composition. This finer grldwork is used to Integrate the non-linear 
ordinary differential equations that describe the chemical reactions that 
occur in the pyrolysis gas. The mass-density equations use the results 
of both the energy equation and chemical kinetics in conjunction with the 
pyrolysis rate to generate Instantaneous density values for all 
components at each coarse grid location. The coarse grid is more than 
adequate for approximating the energy equation. The finer grldwork is 
required to give adequate approximations to the highly non-linear 
chemical-kinetics equations, rhe accuracy to be expected using this 
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numerical approach is indicated by the following examples* 


Comparisons with Exact Solutions 

Two problems for which exact solutions exist were modeled with the 
ablation program and solved numerically. The resulting numerical 
solutions were then compared with the exact solutions. The first problem 
selected was to find the temperature response of a heat sink, with 
constant material properties-, that had a constant-heating rate applied to 
one face and was perfectly insulated on the other face. A comparison of 
numerical and exact solutions to this problem is shown in figure 3. The 
percent error in the numerical solution is plotted for three points 
within the solid: near the surface, the mid-point of the slab, and the 

back surface of the slab. Since the exact solution has a discontinuity 
at time equal zero, the comparison is started at a non-dimensional time 
(kt/l/pCp) of 0.025. This figure shows that 'r^ars in the numerical 

solution are negligible. 


The second problem was to find the temperature response of an 
ablation material undergoing quasi-steady ablation when pyrolysis 
reactions- in-depth are present, and the heat of pyrolysis is zero. The 
zero heat of pyrolysis decouples the mass and energy equations thus 
permitting exact solutions (ref. 8). The resulting exact temperature 
equation is an exponential function of distance. This problem was 
modeled on the computer by imposing constant environmental conditions 
and allowing the quasi- steady conditions to develop as they would in a 
semi- inf inite slab. After about 140 seconds, the numerical solution 

approached a quasi- steady condition by oscillating about the exact 
solution. In the exact solution, the surface recedes at a constant 
velocity, and the temperature and density profiles do not change with 
time. The transient numerical solution oscillated slowly about the exact 
solution. Figure 4 shows the temperature profile, obtained numerically 
after 140 seconds at constant conditions, compared with the exact 
solution. Agreement between the solutions is good. Similar agreement is 
obtained at all greater times investigated. 


COMPUTER PROGRAM 
Using the Program 

To assist the user, three tables are included in this paper. Table 
I contains two sections. The first section is in alphabetical order 
according to physical symbols. These symbols were used in the 
mathematical modeling of the physical system. Each physical symbol is 
defined in section one and the equivalent FORTRAN symbol is given. A 
number is also assigned to each symbol for use as cross-reference to the 


other tables. The second section of Table I contains the same 
information about items that are important enough to require definition 
but were not assigned a physical symbol. 

Table I should be used to go from the mathematical formulation in 
reference 1 to symbols used in this computer program. This table can 
obviously be used as an updated symbol list. 

Table II contains a list of input and output variables in the order 
of appearance in the I/O lists of the computer program. This table can 
be used as a reference when preparing inputs. Two columns of numbers are 
included in this table. The first column, in parenthesis, is the item 
number of that FORTRAN variable. The third column is the number assigned 
to the FORTRAN symbol in Table I. 

The primary method of input is by use of the namelist entitled 
NAME1. The remainder of the inputs are made by individual FORMAT. The 
appropriate FORMAT to use for each input is listed under comments. The 
latter inputs, not included in NAME1, are those that will be changed 
infrequently when working with a particular material system. 

Table III is an alphabetical listing by FORTRAN symbol. The 
majority of these symbols are used internally and have no physical 
significance. When appropriate, the corresponding symbols are listed 
opposite the FORTRAN symbol. The location of the FORTRAN symbols in 
Tables I and II are listed for cross-reference purposes. Table III will 
be of use when working internally with the program. Internal items will 
only be known by FORTRAN symbol. 

A brief description of the table look-up subroutines used in the 
program is included as Appendix A. The inputs for these sub-routines are 
input in standard Namelist form in NAMEl. 

Experience has shown that computing time is virtually unaffected by 
the number of finite difference stations used as long as the total 
(I+J+M) is 100 or below. Consequently, it is recommended that I+J be 
greater than 50 and M greater than 10. It does not save computing 
time to use a smaller number of stations but it will increase truncation 
errors . 


Example Case 

The example case included herein is moderately complex so that a 
number of the more commonly used factors of the program are demonstrated. 
A copy of the actual output is printed in Appendix B. The ablation 
material in this case contains seven solid pyrolysis reactions, all of 
which leave a solid residue upon pyrolyzing. The ablation material is 
modeled by 51 finite-difference stations spaced at 1/n-l or 1/50 • 0.02 
distance apart. The computer automatically calculates the number of fine 
grid points used in integrating the chemical reaction equations. The 
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material is subjected to an environment that is typical of entry at lunar 
return velocity. 

The first section of printout in Appendix B is a list of the 
NAMELIST entitled NAME1. This list continues for 14 1/2 pages and is 
terminated by $END. The values input in this section are typical for 
epoxy ablation materials and can be used where specific data are not 
available. Immediately following NAME1 is a list of the inputs required 
by the pyrolysis-gas kinetics sub-program. The reading and printing for 
this section are done in subroutine IN0UT. The values listed in this 
section are for the chemical species and kinetic equations listed in the 
subroutine. The last section of input is a tabular listing of selected 
data from NAME1. These tabulated items are frequently referred to and 
are relisted for convenience. 

The results of calculations are headed by the title "Beginning of 
Output." The output is keyed to the current computed value of time 
designated as TAU. The first cycle of output occurs after only one time 
step or DELTAU. The first print., in this case, occurs after only 
.0078125 second. As would be expected after such a brief time interval, 
the temperature of the ablator, T, is virtually unchanged from the 
initial condition. Even the surface temperature has risen only 3.84 
degrees. All points of the ablator are at such a low temperature that no 
pyrolysis is occurring. Consequently, the gas flow past all points 
within the solid, WWPDOT, is zero and the density of all reacting 
components of the solid, RHO, remain at initial values. 

One more print sequence is shown at a time of 80.023437. The 
surface temperature has increased to 5400.83 R. Pyrolysis has proceeded 
to a point that the density of all seven solid components has been 
reduced to zero at the surface. Deposition of carbon has proceeded to a 
point that the density at the surface is 23.5441 pounds per cubic foot. 
The average molecular weight of the pyrolysis gases as they enter the 
boundary layer is 11.114. This value is much lower than the average 
molecular weight of 37.815 as produced by pyrolysis. The temperatures 
are so high that a large percentage of the pyrolysis gas should be 
converted to solid carbon. An examination of the rates of deposition of 
carbon confirms this assumption. The sum of the rate of deposition, 

WDEP, of carbon represents 44% of the mass formation rate of pyrolysis 
gas. 


CONCLUDING REMARKS 

The computer program which has been described (CHAP IK) has the 
capability of computing the response of a material to a prescribed heat 
input at the surface. This program is a modification to an existing 
program (CHAP I) . This manual contains a brief explanation of the 
modification that was made to the original program. Use of CHAP III 
should not be attempted without familiarity with the mathematical 
formulation of CHAP I. The tables in the appendix are complete and 


contain all FORTRAN symbols used in the program. A briel explanation 
included for each symbol used in input or output. 


APPENDIX A 



LIBRARY SUBROUTINES 


The following are Instructions on the use of library TABLE LOOKUP 
subroutines contained In the computer program, ^ sheets are cobles 
of pages from the computer programing manual used at the Langley 

Center. 
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SUBROUTINE FTLUP 

NOTE: This subroutine is a special case of MTLUP (El. 5) ond is 

included to provide continuity of existing programs. New 
programs should use MTLUP instead of FTLUP. 

LANGUAGE : FORTRAN 

Computes y = F(x) from a table of values using first or second 
order interpolation. An option to give y a constant value for 
any x is also provided. 

CALL FTLUP (X, Y, M, N, VARI, VARD) 

X - The name of the independent variable x. 

Y - The name of the dependent variable, y = F(x) 

M- The order of interpolation (an integer) 

M = 0 for y a constant as explained in the NOTE below. 

M = 1 or 2, First or second order if VARI is 
strictly increasing (not equal) 

M - '1 or -2. First or second order if VARI is 
strictly decreasing (not equal) 

N - The number of points in the table (an integer). 

VARI - The name of a one-dimensional array which contains 
the N values of the independent variable. 

VARD - The name of a one -dimensional array which contains 
the N values of the dependent variable. 

NOTE: VARD(I) corresponds to VARI(I) for I = 1, 2, ...» N. 

For M = 0 or N ^ 1, y = F(VARI(1)) for any value of x. 

The program extrapolates. 

RESTRICTIONS: All the numbers must be floating point. The values of the 
independent variable x in the table must be strictly 
increasing or strictly decreasing. 


PURPOSE: 


USE: 
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SUBROUTINE DISCOT 


LANGUAGE: 

PURPOSE: 


USE: 


FORTRAN 

SINGLE OR DOUBLE INTERPOLATION SUBROUTINE FOR 
CONTINUOUS OR DISCONTINUOUS FUNCTIONS 

Given some function with two independent variables, 
x and z, this subroutine performs K x th and K x th 
order interpolation to calculate the dependent 
variable. In this subroutine all single line 
functions are read in as two separate arrays 
and all multi-line functions are read is as 
three separate arrays, i.e. 

i - 1) 2, ... L 
Y j . . .M 

Zj^ k = 1 , 2 , . . . N 

CALL DISCOT (XA, ZA, TABX , TABY , TA3Z, NC, NY, 

NZ, ANS) 

XA - The X argument. 

ZA - The Z argument (may be the same name as X on 
single lines). 


TABX - A one-dimensional array of X's. 


TABY - A one-dimensional array of Y's. 
TABZ - A one-dimensional array of Z*s. 


NC - A control word that consists of a sign and three 
digits. The control word is formed as follows: 

(1) If NX = NY, the sign is -. If NX / NY, NX 
is computed by DISCOT as NX = NY/NZ. The 
sign is + and may be omitted if desired. 

(2) A 1 in the hundreds position of the word 
indicates that no extrapolation occurs 
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DISCOT 

above Z max. Kith a zero in this position 
extrapolation occurs when Z >Z max. The 
zero may be omitted if desired. 

(3) 1-7 in the tens position of the word indi- 
cates the order of interpolation in the 
X direction. 

(t) 1-7 in the units position of the word 

indicates the order of interpolation in 
the Z direction. 

NY - The number of points in the Y array 

NZ - The number of points in the Z array. 

ANS - The dependent variable Y. 

The following programs will illustrate various 
ways to use DISCOT. 

Case I. Given Y = f(x) 

NY = 50 

NX ( number of points in X array) = NY 
Extrapolation when Z > Z max 

Second order interpolation in X direction 
No interpolation in Z direction 
Control word = -020 

1. DIMENSION TABX (50), TABY (50) 

1 FORMAT ( 8e 9.5) 

READ (5, 1) TABX, TABY 
READ (5, l) XA 

CALL DISCOT (XA, XA, TABX, TABY, TABY, 
-020, 50, 0, ANS) 

Case II Given Y = f(x, z) 

NY = 600 
NZ = 10 

NX = NY/NZ (computed by DISCOT) 
Extrapolation when Z >Z max 
Linear interpolation in X direction 
Linear interpolation in Z direction 
Control word = 11 
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DIMENSION TABX (800), TABY (800), TABZ (10) 
1 FORMAT ( 8E 9.5) 

READ (5, 1) TABX, TABY, TABZ 
READ (5, 1) XA, ZA 

CALL DISCOT ( XA , ZA, TABX, TABY, TABZ, 

11, 800 , 10 t ANS ) 


Case III. Given Y = f(x, z) 

NY = 8 00 
- NZ '■* 10 
NX = NY 

Extrapolation when Z >Z max 

Seventh order interpolation in X direction 
Third order interpolation in Z direction 
Control word = -73 

DIMENSION TABX ( 800 ) , TABY { 800 ), TABZ ( 10 ) 

1 FORMAT ( 8E 9.5) 

READ (5, 1) TABX, TABY, TABZ 
READ (5, 1) XA, ZA 

CALL DISCOT ( XA , ZA, TABX, TABY, TABZ, 

-73, 800, 10, ANS) 

Case IV. Same as Case III with no extrapolation above 
Z max. Control word = -173 


CALL DISCOT (XA, ZA , TABX, TABY, TABZ, -173, 
800, 10, ANS) 


RESTRICTIONS: See be of METHOD for restrictions on t abulat i ng 

arrays and discontinuous functions. The order 
of interpolation in the X and Z directions may 
be from 1-7. 


METHOD: 


The following subprograms are used by DISCOT: 
UNS, DISSER, LAGRAN. 


Lagrange’s interpolation formula is used in both 
the X and Z direction for interpolation. This 
method is explained in detail in Methods in Numeri- 
cal Analysis by Nielsen. The search in both the 
X and Z direction observe the following rules: 


1. X < Xp the routine chooses the following points 
for extrapolation. 
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I 
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$ 

t 

? 

i 

r 

! 


f 

* 


^1* X 2* • • * • Y i» Y 2* * * * * Y k + 1 

?, X N X the routine chooses the following 

n 

points for extrapolation. 

x n-k> X n-k+l» * * * * X n Y n-k» Y n-k+i* * * ln 

3. X 1 X n the routine chooses the following 
points for interpolation. 

k is odd x i _ k+1 . X i. k+ i+l> * * * X i-k+l+k 

2 2 2 

Y j. .k+i » Y i-k+l_+l 5 * ' * * f i-k + l+k 

2 2 2 


\ 

{ 

i 

\ 


i 


i 

! 


k is even X i-k , X i _ k+] _ J x i-k+k 

2 2 2 


Y i_ k * Y i-k+l’ 1 i-k + k 

2 2 2 


4. 


If any of the subscripts in Rule 3 become 
negative or greater than n (numcer of points 
Rules 1 and 2 apply. When discontinuous 
functions are tabulated, the independent 
variable at the point of discontinuity is 
repeated, i.e. 


k - 


2 (X 1$ Xg ,X^ ,X^ ,X^ , ,Y 2 ,^3 » ^ ,Yr; , Y, 


The subroutine will automatically examine th< 
points selected before interpolation arid if 
there is a discontinuity, the following rule 
apply. Let X d and X d+1 be the point of 
discontinuity. 
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PROGRAM ING 
MANUAL 




a* X < points previously chosen are modified, 
for interpolation as shown 


X d-k> X d-k+l» * * * * X d 
Y d-k » Y d-k+l* • • * • Y d 

b. X >X d points previously chosen are modified for 
interpolation as shown 

X d+l J X d+2 * * * * ‘ X d+k 

Y d+1» Y d+2* * * * * Y d+k 

c. When tabulating discontinuous functions, 
there must always be k+1 points above and 
below the discontinuity in order to get 
proper interpolation. 

When tabulating arrays for this subroutine, 
both independent variables must be in 
ascending order. 

In some engineering programs with many tables, 
it is quite desirable to read in one array of 
x's that coula be used for ail lines of a 
multi-line function or different functions. 

The above not only saves much time in preparing 
tabular data, but can also save many locations 
previously used when every y coordinate had 
to have a corresponding x coordinate. Even 
though the above is not always applicable, 
the subroutine has been written to handle this 
si tuation. 



Another additional feature that may be useful 
is the possibility of a multi-line function 
with no extrapolation above the top line. 
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ACCURACY: 

REFERENCE: 

STORAGE: 

SUBPROGRAMS 

USED: 


OTHER CODING 
INFORMATION 

SOURCE: 

RESPONSIBLE 
PERSON : 


A function of the order of interpolation used. 
Nielsen, K.L.; Methods in Numerical Analysis 
DISCOT - 555g locations 


UNS 1 0 locations 
DISSER 110§ locations 
LAGRAN 55g locations 


NONE 

SHARE Library, General Motors Corp» , Allison Div, 


Vivian P. Adair 
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***** EXAMPLE CASE 



8KAME1 


2 

TZ 

HCP 

NCP 

TtABCP 

CPTAB 

PCPP 

NCPP 

TT8CPP 

CPPTAft 


0. 54E *03 1 

0.54E*03, 

0.546*03, 

0.546*03, 

0.54E*C3, 

0.54E*03, 

0.546*03, 

0.54E*03, 

0.54E*C3, 

0.546*03, 

0.546*03, 

0. 54E*03, 

0.54E*C3t 

0. 54E*03 , 

0.54E*C3, 

0.546*03, 

0.546*03, 

0. 54E*03, 

0.546*03, 

0.546*03, 

0.54E*C3, 

0.546*03, 

0. 54E *C3, 

O. 54E*03» 

0. 546*03 , 

0.54E*03, 

0.546*03, 

0.546*03, 

0.546*03, 

0.546*03, 

0.546*C3, 

0 . 54E *03, 

0. 546 ♦ C3 , 

0.546*03, 

0.546*03, 

0.546*03, 

0.546*03, 

0.546*03, 

0.54E*C3, 

0.546*03, 

0. 54E*C3, 

0.546*03, 

0.546*03, 

0.546*03, 

0.54E*C3, 

0.546*03, 

0.546*03, 

0.54E*03, 

0.546*03, 

0. 54E *03 , 

0.54E* *.3, 

0.546*03, 

0. 546*C3, 

0.546*03, 

0, 54E *C3 , 

0.54E*03, 

0 .546*03, 

0.546*03, 

0.546*03, 

0.54E*03, 

0. 546 *C3 , 

0.0, 0.0, 

0.0, 0. 

0, 0.0, 0. 

0 , 0.0, 0 

. 0 , 0.0, 

0.31686*00 

0* Of 0* 0 1 

t C.3021E*00, C.277lE*00, 0.25216*00, 0. 

C.Ot 0.0, O.C, 

16256*00, 

0*0i C* 0 1 

0.0, 0.0 

t 0.0, 0 

• Ot 0*0 • 0 

•Ot 0« Of 

0.0, 


■ i» 


» 8 1 

_ 0.726*03, 0 • IC8E*04» 0.1446*04, 0.186*04, 0.2166*04, 0.2526*04, 

0.25746*04, C.7E+04, O.Ot O.Ot O.O* 0.0, 0.0, O.Ot O.Ot 


0.0, 0.0, 

0 « Cf 0 » 0 9 0 • C 9 

O.Ot 

0.0, 

0.0. 

O.C, 

0.0, 


• 0. 25C4E*0l), 

C#299iE ♦■00» C. 

3478E*00 

1, o 

• 3965E *00, 0. 

4452E*00, 

0. 49396*00, 

0#5012E*00, C. 

5E*00t 

0.0, 

0.0, 

0. 0, 

0.0, 

0.0, 

0.0, 0.0, 

0*0, 0*0, 0*0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 


O.Ot 
» 0 , 

■ 0 » 

• O.Ot 0. 56E*03» 0.666*03, 0.766*03, O.B66*C3, 0.96E*03, 

0.1066*04, 0. 1166*04, 0.121E*04» O.Ot O.Ot O.Ot O.Ot O.Ot 

O.Ot O.Ot O.Ot O.Ot O.Ot O.Ot O.Ot O.Ot O.Ot O.Ot O.Ot 


• 0.36*03, 0. 3296*00, 0.364E*00, 

0.4246*00, 0.425E*00, 0.4256*00, 


0.39 7L*0C, 0.4066*00, 0.4186*00, 

0.0, O.Ot O.Ot O.Ot 0.0, 




0.0, 

0.0, 

0.0, < 

MCPOP 


0, 



NCPOP 

» 

0, 



fTCPDP 

U 

0.0, 

0.0* 

O.C* I 



0.0, 

0.0* 

0.0, I 



0.0* 

0.0* 

0.0, 

cpopie 

m 

0. 116*30 , 

O.C, 0 



0.0, 

0.0* 

O.C* 



0.0, 

0.0* 

0.0* 

PCPBAR 

m 

1, 



NCPBAR 

m 

2, 



TTBCPB 

m 

0. 1E*0 A, 0 

.26*04, 



0.0, 

0.0, 

0.0, 



0.0, 

0.0, 

c.c. 

cpbtab 

m 

0. IE *04, 0 

.26*04, 



0.0, 

C.O, 

0.0. 



0.0, 

0.0, 

0.0, 

PK 

m 

1, 



AR 

m 

8 * 




O.Ct O.O* 0.0* 0.0* 

0 • C * 0.0 1 0*0) 0.0* 


0 . 0 * 0 . 0 * 0 . 0 * 

0 . 0 * 0 . 0 * 0 . 0 * 


0 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 0 . 0 , 


0 . 0 * 


0 . 0 * 0 . 0 * 0 . 0 , 0 . 0 ) 


* 0 . 0 , 0 . 0 , 0 . 0 , 

0 . 0 * 0 . 0 * 0 . 0 * 


TTA8K * 0.5E*03» 0.1E*04, 0.15E*04, 0.26*04, C.25E*04, 0.3E*04, 

0.35E*04, 0.46*04* 0.0* 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

0.0, C. 0, C.O, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0* 


CKTAB • O.14E-04, 0.14E-04, 0.388E-04, 0.813E-C4, 0.117E-03* 0.167E-03* 

0.2E-03* U.2E-03* 0.0* 0.0* 0.0, 0.0, 0.0. 0.0* 0.0* 0.0* 
0.0* 0.0, 0.0* 0.0* O.C* 0.0, 0.0, 0.0, 0.0, 

MCKNP * 0, 

NCRNP ■ 8* 

TTCKNP * 0.5£*C3* 0.1E*04, 0.156*04, 0.2E*04, 0.25E*C4 t 0.3E*04, 







0.35E»C4, 
O* Of 0# Of 

C.4E+04, 0. 

0.0, 0. 0, 

Of 0.0* O.Of Of Of Of Of Of Of OfOf 0.0* 
Of Of 0* Of 0# 0 1 0«0f O.Of 

CKNPTB 

• 0.14E 

0.2E- 
0.0, 

-04, 0.14E-04, 0 

03, 0.0, 0.0, 0 

0. o, o.o, 0.0, 

.3888-04* 0.8138—04* 0# 1 i 7E-03# 0.1678-03* 

•Of O.Of O.Of O.Of O.Of O.Of O.Of O.Of 
O.O* 0. Of 0.0* O.Of 

HCKNOP 

■ 0, 




NCKNDP 

* 0, 




TCKNOP 

* 0.0, 

0.0, 

0.0, 0.0., 

0.0* 

CKNDPT 

• 0.417E-04, 

o.o« o.o, 
0.0, 0.0, 
0» C> 0.0, 

0 u 0, 0. 0 , 

0.0, 0.0, 

0.0, 0.0, 

0.0, 0.0, 

0.0, 0.0, 

0.0, 0.0, 

O.Cf O.Of O.Of O.Of O.Of O.Of O.Of 

O.Of O.Of O.Of O.Of O.Cf O.Of O.Of 

O.O* O.Of O.Of O.Of OfOf O.Of O.Of 

O.Of O.Of O.Of O.Of O.Cf 0.0*. O.Of 

O.Cf O.Of 0.0* 

HHW 

% 1, 




NhU 

• 14, 




TfABHW 

* 0.5£»03, 

0.35EMJ4, 
0.65E*04, 
0.0, 0.0, 

0.lE*-04, 0. 15E»C4, 0.2E+04, C.25E*04, 0.3E«>04, 

0.4E*04» 0.45E *-04, 0.5E*04, 0.55£»04, 0,6E*04, 

0.7E*04, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

0.0, 

HWTAB 

» -0. 1E*02, 
0.818E+03. 
0.206E >04, 
0.0, 0.0, 

0. 111E*03, 0 

0.98E«-03, 
0.255E*04, 

0 • 0 9 0 • 0 f 

i.239Et03, 0.375E*C3, 0.517E*03, 0.665E*03, 

0,117 3E*04» 0.13B4E*04, 0.17E*04, 

0.3E»04, 0,0, 0,0, 0.0, 0*0, 0*0, 0.0, 

0.0, 

MALPHA 

* 0, 




NALPHA 

■ 0, 




TALPHA 

* 0.0, 

0.0, 

0 • Of 0*0 1 

o.o. 

ALPHAT 

* 0. 1159E+01, 0.0, 0.0, 

0,0, 0.0, 

MALPH1 

■ 0, 





ir- 



( 


s 








NALPH1 - 0, 

TTALl - 0.0, 0.0, 0.0, 0.0, 0.0, 

AUTAB « 0.5E*00, 0.0, 0.0, 0.0, 0.0, 

HALPH2 >0, 

NALPH2 * 0, 

TTAL2 - 0.0, 0.0, 0.0, 0.0, 0.0, 

AL2TAB - 0.1E+01, O.C, 0.0, 0.0, 0.0, 

MhE • 0, 

KHE * 0, 

TTABHE - 0.0, 0.29E*02, 0.39E*02, C.49E*02, 0.59E*02» 0.79E*02» 

0.59E*02» 0. 1 19E*03, 0.139E*03, 0.149E*03, 0.159E*03, 


0. 179E *03 , 

0.199E*03, 

C.219E*03, 

0. 239E*C3 , 

0. 259E *03 , 


0.295E *03, 

0. 329E *03 , 

0.3696*03, 

0.399E*03, 

0. 419E *03 , 


0.445E+03, 

0.489E* 03 , 

0. 505E*03 

0. 529E *03 , 

0.5496*03, 


0.565E*03, 

0 ‘>996*03 , 

0.11E*04, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 0.0, 

O.u, 0.0, 

0.0, 0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 0.0, 

0.0# 0.0# 

O.C# 0.0# 

0.0, 

0.0, 

O.C, 

0.0, 

0.0, 

0* G» 0.0# 

O.C# 0.0# 

0. 0 , 0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

O.C, 0.0, 

0.0, 0.0, 

0.0# 0. 0# 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

O.C# 0. 0 # 

0.0# 0.0# 

0.0, D. 0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

O.C# 0.0# 

0.0# 0.0# 

O.C# 0.0# 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 



HETA8 * 0.19226*05. C.L922E*05, C.1922E*05, 0.1534E*05, 0.1934E*05, 

0.1922E*05, C. 18736*05, C.1764E*05, 0.1565E*05, 0.13186*05, 

0.11046*05, C.94E*C4, 0.8C8E*04, 0.7046E*04, 0.61886*04, 

0.54781104, 0.488EO4, 0.4375E*04, 0.3948E*04, 0.3575E*04, 

0.3256E*C4, C.2977E04, C.2724E*04, 0.2491E*04, 0.2273E*04, 

0.2CCEO4, G. 18 55E 04, 0.1649EO4, 0.14386*04, 0.1232EO4, 

0.10256*04, C.832E + 03, 0.65676*03, 0.504E+03, 0.3767E*03, 

0.27cB i :t03» C.2112E *03, C.1404E*03, 0.1404E03, 0,9768E*02, 

0.57716*02, C.4682 E *02 , 0.3277E*02, 0.229E*02, 0.162E*02, 

»s» 0.12CE*02, 0.586 *01 , 0.72E*01, 0.5E*01, 0.8E-02, 0.8E-02, 






o £ 
£ 


§6 



TTABQC 


QCTA8 


0.6E-02, 0.86-02, 0.86-02, 0.8E-02* 0.8E-02* 

G.8E-02, 0.8E-02» 0.8E-C2, 0.8E-02, 0.8E-02, 0.8E-02, 


0. Ot 0. 0 , 

0 . 0 . 0 . 0 . 

0 . 0 . 0 . 0 » 

0.0. 0.0. 


0 . 0 . 0 . 0 . o.c, 
0.0, O.C, O.c, 
0.0, 0.0, 0.0, 


1 0 , 0. 
.76*02, 
,136*03 
. 19E *-03 
,256*03 
.316*03 
> 376*03 
.436*03 
.496*03 
.556*03 
.616*03 
. 0, 0. 
.0, o. 

. C , 0 . 


1E+02 , 0.2E*02» C.3E*02» 0.4E*02» 0.56*02, 0.6E*02» 

0 . 8E*02 , 0.9E*02, 0. IE*03» 0.11E*03, 0.12E*03» 


0.146*03, 
0.2E*03. 0 

0.266*03 , 

0 . 3 2E *-03 , 
0.38E*03, 

0 • 446*03 , 
0.56*03, 0 

0. 566*03 , 

0 . 1 3E *04 , 
0 . 0 , 0 . 0 , 
0 . 0 , 0 . 0 , 
0 . 0 , 0 . 0 , 


0.156*03, 0 

0.216*03, 0. 

0.276*03. 0 

0. 33E*03» 0 

0. 39E*03, 0 

0.45E*03, 0 

0. 51E*C3» 0. 

0.576*03, 0 

0.2E*C4, 0. 

1, 0.0, 0.0, 

i, O.C, 0.0, 

l , 0*0, 0.0, 


0.166*03 , 0.176*03, 0.18E*03, 

C.22E*03, 0»23E*03» 0.246*03, 

0.28E*03 , 0.29E *03, 0.3E*03» 

0. 34E*03 , 0.356*03, 0.36E*03, 

0.46 *03 , 0.416*03, 0.426*03, 

0.466*03 , 0.47E*03, 0.486*03, 

0.526*03, 0. 53E*03, 0.546*03, 

0.586*03, 0. 59E *03, 0.6E*03, 

0.0, 0.0, 0.0, C.O, 0.0, 0.0, 

0, 0.0, 0.0, 0.0, 0.0, 0.0, 

0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

0 , 0 . 0 , 0 . 0 , 


0. 365E *01 , 
0.1927E*03, 
0, 33376*03, 
0.9796 *02 , 
0.472E+02, 
0.2*56*0 2, 
0.1*46*02, 
0.22L*0l, 
0.iE*00» C 
0.5E-01, C 
0. 26-01, C 
0. IE-02 , C 


0.756*01, 0. 168E*02, 0.406E*02, 0.9616*02, 

, 0.31776*03, 0.4385E*03, 0.49016*03, 0.43446*03, 

, 0.26186*03, 0.16986*03, 0.13976*03, 0.11636*03, 
0.8316*02, 0,7126*02, 0.615E*02, 0.533E*02, 

0.422E*02, 0.3826*02, 0.^496*02, 0.3216*02, 

0.2716*02, 0.24£E*02, 0.2186*02, 0. 1876*02, 

0.1226*02, 0.926*01, 0.67E*C1, 0.476*01, 0.336*01, 

0.226*01, 0. 1 46 *C1 , 0.1E*01, 0.6E*00, 0.46*00, 

0.26*00, 0.156 *CC , 0.1E*00, O.BE-Ol, 0.66-01, 

0.4E-01, 0.36-01, 0.3E-01, 0.36-01, 0.2E-01, 

0.2E-01, 0.2E-01, 0.2E-01, 0.16-01, O.IE-Ot, 0. 16-02, 

0.1E-02, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

0.0, 0.0, O.C, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

n . a. 0 . 0 . 0.0, 0.0, 0.0, O.C, O.C, 0.0, 0.0, 


0.0, 

0.0, 

0.0, 

0.0, 

o.c, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

o.c, 

0# 0 1 

0.0, 

o.o, 

0.0, 

0.0, 

o.c, 

O.Ot 


* 0 , 




KinvnS w 
si aova rrvi 


* 0 


© © 

«g 


TTABQR ■ O.O* 0.3IE+02, 0.36E+02, C.41E+02, 0.46E+C2, 0.51E+02, 

0 • 56£ *- 02 1 0.59E+02, 0.61E+02. 0.63E+02, 0.65E+02, 0.79E+02, 

0.81E+02, 0 • 1 it +02 , 0.91E+C2, 0.261E+03, 0.291E+03, 0.*41E+< 

0.461E+03, 0.491E+03, 0.541E+03, 0.LIE+04, O.O, O.O, O.O* 

O.Of 0.0. 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, o.o, 


0.491E+03, 

0. 54 LE+03, 

0.LIE+04, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 
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ACOEF 


O.C, 

0.0, 

0.0, 










VOEL 


0 . 0 , 



ELAMTS * 0.736400, 0.756400, 0.75E400. 0.75E*00» C.75E400. 0.75E400, 

0.75E*-0C» 0.75E ♦00» 0.75E4C0. 0.75E400. 0.75E4C0, 0.75E400, 

0. 756*00, 0.75E400. 0.75E4C0. 0.75E400. 0.756400, 0.75E*00. 

0.75E40C. 0. 75t*-00* 0.75E*C0, 0.75E400, 0.75E*00» 0.75E400, 

0. 75t *-00 * 0.756*00, 0.75E4C0, 0.75E*00» 0.75E4C0, 0.756400. 

0.75E*OC. 0.15fc»00* 0. 756*00. 0.75E*00, 0.756400. 0.75E*00, 

0.75E4CC, 0.7>C*-00» 0.75E*C0, 0.75E400. 

TTELAM » 0.$E*03, 0.56*03, 0.5E*03, 0.5 p -03, 0.56*03, 0.56*03, 0.56*03, 

0.SE*03, 0.5e*03» 0.1E*C5» 

PEL AM » 0.1E-01, 0.5E-01, O.lE*01, 0.26*01, 

IPRIKT * 0. 

KT.XDP » 0# 

ICP6AR • 2* 

CPPT7 » 0* 3E *-00 , 0.0. 0.0. 0.0. 0.0. 0.0. 0.0. 0.0* 0.0. 0.0. 

TTHCCMB » 0.18E*04, 0.276*04, 0.36£*C4, 0.456*04. 0.54E*04, 0.636*04. 

0.726*04, 

HCCM6TB ■ 0.41 1E*C4, 0.42666*04, 0.44546*04, 0.48716*04, 0.6265E*04, 

0. 1022E*05» C.1354E*05» 0.4liE*C4, 0.4266E*04, 0.44476*04, 
0.46976*04. 0.52956*04, C.6995E*04, 0.13C5E*05» 0.411E*04, 

0.4266E*04, 0.44466404. 0.4656E*04, 0.49836*04, 0.56796*04, 

0. 71346*04, 0.4116404. 0.4266E*04, 0.4445E404, 0.4643E404. 
0.48646*04, 0.5245E404, C. 98696404. 

PHCOMB * 0. IE 40C . 0. IE 401 . 0 . IE 402 * 0.1E403. 

MHCCM9 « $6862409532376811. 

NPHCCMB ■ $6862409532376634. 

MCP7 • 0, 

NCP7 * Or 

T7CPP * O.C. 0.0. O.C. 0.0. 0.0. 0.0. 0.0. 0.0. O.C, 0.0. 


CMASSF 

CHANGE 

DOELHP 

CSOL 10 

RHOZ 

NPVR 

AP 

BP 

RHO 


0.52£*C0, 

C.2E*01, 

O* 125E*03, 0.1256*03, 0. 1256*03, 0.166E*03, 0.166E*03, 

0.166E,03, 0.276*03, 

0* 25E*CC , 0.25E*C0, 0.256*00, 0.25E*00, 0.256*00, 0.256*00, 

0.256*00, 


0.60066*00, 

C.7856E *00, 

0.101666*02, 

0.800256*00 

, 0* 136125E< 

0.608856*01, 

0.1326*02, 




7, 

0.51666666*09 

, 0,256*06, 

0.216666676*08, 0. 2 166 6667E*06, 

0. 966666E *07 , 

0. 136*11 , 

0.83333336*15 

* 


0.6967 185E*05 

, 0.6366919E 

*05, 0,602266E*05, 0.3011336*05, 

0. 526831 86*05 

, 0.73992686*05, 0.998060896*05, 


0. 60066*00, 

0.60066*00, 

C. 60066*00, 

0.60066*00, 

0.60066*00, 

C.60C6E*00, 

0 • 6006E *00, 

C.60C6E*00, 

0,60066*00, 

0.6006E *00, 

0. 60C66*00, 

C .6006E *00, 

C . 60Ct>E*00 , 

C . 6006E *00 , 

0. 6006E *00, 

0. 60C6E*00, 

C. 60066*00, 

C.6006E*00, 

0.6006E*00, 

0.60066*00, 

0. 60G66*00« 

C • 6006E *00 , 

C* 60C6E *00 , 

C.6006E*00, 

0.60066*00, 

0. 60C66 *00 , 

C .6006 E *00 , 

C. 60 066*00, 

0.6006E*00, 

0.60066*00, 

0.60066*00, 

0.60066*00, 

C.60C6E*00, 

0.6006E *00 , 

0.60066*00, 

C. 60C6E *00, 

0 .60066 *00 , 

C.60C6E*00, 

0. 60Q6E *00, 

0.60066*00, 

0, 6QC6E*G0» 

C. 60066*00, 

C.6006E*00, 

0,6006E*00, 

0.60066*00, 

0.60C6E*00, 

C .60066 * 00 , 

C.fc0C6E*00, 

0, 60C6E *00, 

0,6006E*00. 

C. 60C.6E *00, 

C .6006 E * 00 , 

C.60G6E *00, 

0, 6006 E *00 , 

0.6006E*00, 

0. 60C6E*00, 

0.60066*00, 

0.60C6E*00, 

0,60066 *00, 

0.6006E *00 , 

C.60C6E*00, 

0.7856E *00, 

C. 78566*00, 

0.78566*00, 

0, 7856E*00» 

0.78566*00, 

0.7856E*00, 

0.78566*00, 

0.7856E*00, 

0.78566*00, 

0. 7856E*00, 

0.7656E *00, 

C.7856E*00, 

0. 7856E *00, 

0. 7856E*00, 

0. 78566*00, 

C.7856E*00» 

C.7856E*00, 

0. 7856E *00, 

0.78566*00, 

0. 78566*00, 

0 .7856E * 00 , 

C.7856E*00, 

0.78566*00, 

0.78566*00, 

0. 7tt56E*C0, 

0.7856F*00, 

C. 7856E *C 0, 

0,7 856E *00, 

0. 7856E*00, 

0. 78566*00, 

0,7356 fc*00. 

0.7856E*00, 

0.78566*00, 

0. 7tt56E*00, 

0. 78566*00, 

0.7 356E *00, 

0. 78566*00, 

0.7856E*C0, 

0.78566*00, 


(ORIGINAL PAGE IB 
OF POOR QUALITY 


0.7854E600, 
0* 7854E* 00* 
0.7854E600, 
0.7854E6C0. 
0. 7854E * 00, 
0. 1 0 16 4E *-02 , 
0. 10 164E *-02 , 
0. 10164E6C2, 
0. 1 C16 £ *02 i 

o. ioio^e *02, 

0. 1G164E*Q2, 
0. 10164E602, 
0. 10 16 4E *02 * 
0. 10164E6G2* 
0. 101646*02, 
0. 10164E*02, 
0. I0164E *02, 
0* tO 02 5E *• 00 » 
0.80C25E*00, 
0.fc0C25E*00, 
O.EOC25E*UO, 
0. 8002 5F *00 » 
0* C0C2 5E *-00 f 
0.60025E*00, 
C.80C25F*OC, 
0. 600256*00, 
0. 8002 * E *-00 * 
C.eOC25E*0O, 
0.80C25E*0G. 
0. 1361256*01 
0. 1.3tl25E*0l 
0. I36125E*01 
0. L361 2 5b *0 1 
0. 136 12 5c *0 1 
0. 136 1 25E *0 1 
0. 136125E*01 
0.136125E*01 
0. l3t!25E*01 
0. 136 l 25E *0 1 
0. 1361256*01 
0. t/t 1 2 5 c ♦ 0 1 
0. 13c 125E »C1 


0.7854E*00, 0. 

C.7854E*00, C. 

0.7B54E600, 0. 

0 • 7 854E *-00 * 0. 

0.78546*00, C. 
G-. 101 64E *02 ♦ 

0. 101 6 4E *02 * 

0. 10164E*32, 

0. IC1 64E *-02 « 

0. 10164E*02, 

C. 101 64E *02 * 

0. 10 164E *-02 • 

0. 10164b *-02, 
0.10164E*-02, 

0. 10 1 64E *-02 r 
C. 10164E*02, 

0. 101646*02, 
C.80025E600, 
C.80025E*OC, 
0.80025E*00, 

C. 8002 5 E *00 1 
C.80025E600, 

0. 800255*00, 

C. 8002 5E *00 « 
0.800255*00, 

0. 8002 5E * 0 C i 
0. 8002 5E *-00 » 

C* 8002 5 E *-0 0 * 
0.80025E*00, 

0. 136 1 25E *01 
0 » 136 1 2 5E *0 1 
0. 136 1 25E * 0 l 
0. 136 1 25E *0 1 
0 • 136 1 25E ♦ 0 l 
0. 136 1 25E *0 1 
0.136125E*01 
0.136125E*01 
0.136125E*01 
0. 136 1 2 5E *-0 1 
0. 136125£*01 
0. 136125E*C1 
0. 1361 25c*01 


7854E*00, 
7854E*00, 
7854E *00, 
7854E*00, 
10164E602. 
0. 10164E*02 
C. 10164E*02 
C. 10 16 4E *-02 
0. 10164E+02 
0. 10L64E602 
0.10164E*02 
0. 10164E*02 
0. 10164E602 
C.10164E602 
0. 10 164E*02 
0. 10164E+02 
C. 1 0 1 6 4 E *02 
C.80025E*00 
C. 80025E *00 
C.B0025E ♦•00 
C.80025E*00 
C.83C25E600 
C.80025E*00 
C. 8002 5E *00 
C.80025E*00 
C.80025E*00 
C.80C25£*00 
C.80025E*00 
C. 600256*00 
0. 136125E 
0. 136125E 
0. 136125E 
0.136125E 
0. 136125E 
0. 136125E 
0. 136125E 
0. 136125E 
0.136L25E 
0. 136125E 
0.136125E 
0. 13612 5E 
0* 136125E 


0.78546*00, 0. 

0*7854E*00, 0. 

0.7854E*C0, 0. 

0 .78546*00 1 O. 
0. 10 164E 602 , 

0. 10164E *02 * 
0. 10164E *02 » 
0. 10164E *02 t 
0. 10164E *02, 
0. 10164E 602, 
0.10164E 602, 
0. 101646*02, 
0.10164E*02, 
0.101646*02, 
0. 10164E *02, 
0. 10164E 602, 
0.80C25E *00, 
C.BC025E 600, 
C. 80C25E *00, 
0.800256*00, 
0.80025E 600, 
0. 80025E 600, 
C.80025E600, 
C* 80025E *00 , 
C, 00025E *00, 
0.80025E 600, 
0.80025E600, 
C.80025E600, 
C • 80025E 630, 
601, 0.136125E 

601, 0. 136125? 

601, 0.136125E 

601, 0 a 136 125E 

601, 0.136125E 

601, 0.136125E 

601, 0.136125E 

601, 0.136125E 

601, Oa 136125E 
601, 0.136125E 

601, 0.136125E 

601, 0.136L25C 

601, 0.13612SE 


7854E600, 

7854E 600 , 
7854E600, 
7854E600, 

0. 10164E602, 
0.10164E602, 
0.10164E602, 
0. 10164E602, 
0.101646*02, 
0. 10164E602, 
0. 10164E 60 2, 
0. 101646*02. 
0.10164E602, 
0.10L64E602, 
0. 10164E602, 
0.10164E 602, 
0.800256*00, 
0. 80025E 600, 
0.80025E600, 
0 • 8002 5E 600 * 
0. 800256 600, 
0.80025E600, 
0* 8002 5E 600, 
0.80025E600, 
Oa 80025E»00, 
0.80025E 600 , 
0. 60025E 600, 
0 . 8002 5E *00, 

* 01 , 

601 , 

601 , 

* 01 , 

* 01 , 

601 , 

* 01 , 

601 , 

601 , 

601 , 

* 01 , 

* 01 , 

♦ 01 , 






o 



RHOC 


0. 1361256*01 , 
0 • 136 125E *0 L . 
0. 136125E*01 . 
0.6088 5E *-01 , 
0.608856*01, 
0.60e85E*01. 
0*60 88 5E *0 1 * 
0.606856*01 • 
0.608856*01, 
0.603856*01 » 
0.6C685E *01* 
0.60t'85E*01 , 
0. 60 88 5E * 0 l , 
0.60 8856*01, 


0.1326*02 
0.1326*02 
0.132E*02 
0.1326*02 
0» 13<1E*02 
0.1326*02 
0. 1326*02 
0. 132E*02 
0. 1326*02 
0. 1326*02 
0. 132E *02 
0. !32E*u2 


0. 1361256 *01 « 
O.X36125E*01t 
0.608856*01, 
0.60885E*CX» 
0.608856*01 , 

0. 6C 8856 *01 * 
0.60865E*01, 

0. 608U5E *01 * 
0.608656 *01 , 
0.60685E *01 , 

C. 608856*01, 

0. 6088 5 E *01* 
0.606856*01, 

0.608656 *01 . 




0. 136125E*01 ♦ 0® 136 125E *0 1 ♦ 

0® 136 125E*0 1 * 0® 136 125E *01 * 

0.60885E*01. O.6O805E*O1, 

0® 60805E *01 • 0.60885E*01, 

0.608856*01, 0.608856*01, 

0.608856 *0 1 1 0.6088 5E*Ol v 

0.60885E *01 , 0.60885E*01® 

0.606356 *01 » 0.608856 *0 it 

C. 6C8 856 *01 » 0.608856 *01 * 

0.60 885E *01 » 0.608856 *0 1 » 

0.608856*01 0® 6088 5E *01* 

Q.60885E *01® 0.60885E*01» 

0.608856*01, 0«60885E*0 1 , 

0® 60885E *01 , 0. 60885E *01, 

0* 132£*02, 0.1 32E*02, 


C.6C885E*01 , 
0.60885E*01, 
C.6C883E*0l, 
0.60885E*01, 
0.608856*01 , 
C.60685E*01, 
0. 60 80 5 E* 0 1 , 
C.6CE85E *01, 
C. 608856*01, 
0. 60885E*01 , 
0.6C885E*0l, 
0.6C885E*0l, 


0 . 0 , 

0 . 0 , 

0 . 0 , 

0.3, 

0.3, 

0 . 0 , 


0.0 

0 



0. 1326*02, 

0. 132E*02, 

0.132£*C2, 


0. 1326*02, 

0 • 132E *02 , 

0. 132E *02 , 


0. 1326*02, 

0.1326*02, 

0. 1326*02 , 


0.1326*02, 

0. 1226*02, 

0. 1326*0 2, 


0. 132E*02 , 

0.1326*02, 

0.1326*02 > 


0, 132-*02, 

0. 132E*02, 

0.1326*02, 


0. 132E*0 2, 

0 . l'ii'E * 02 , 

0. 132E *02, 


0. 132E*02, 

0. 132E*02, 

0. 1326*02 » 


0. 132E*02, 

0. 132E*02. 

0.1326*0 2, 


0 . 1326*02, 

0. 132E*02 , 

0.1326*02, 


0. 132E *02, 

0. 132E*02, 

0. 132E *C.i 


0. 132E*02, 

0.1326*02, 

0.1326*02, 

), 

0 • 0 9 0*0 » 

C • C 9 0 • 0 9 

o.o, 0 . 0 * 

I, 

0 * 0 1 0« 0 1 

0 • 0 1 0 * 0 » 

0 . 0 , 0 . 0 , 

* 

i f 

OoOl 0#0t 

0 • C 9 0# 0 f 

0 . 0 , 0 . 0 , 

), 

0 * 0 9 0» 0 9 

0 • C 9 0 • 0 9 

0 . 0 , 0 . 0 , 

), 

0.0 9 0* 0 9 

0 • C 9 0 .O 9 

0.0, 0.0, 

), 

0 • 0 9 0.0* 

0 » 0 9 0*0 1 



0. 132E *02 , 
0. 1326*02, 
0. 132 E *02 , 
0. 132E *02 » 
0. 132E*02 , 
0. 132E*02 , 
0. 132E *02 , 
0. 1326*02, 
0. 132E *02 , 
0. 132E*02» 
0- 132E*02, 


0.0, 0.0, 0.0, 
0.0, O.O, 0.0, 
O.C, 0 . 0 , 0 . 0 , 

0.0, 0.0, 0.0, 
0.0, 0.0, 0.0, 


ORDER 

* 0.16*01, 0. 

1E*01 

* 0. 

16*01 , 

0. LE*01 , 0 

• 1E*01 

f. 0. 

LE*01, 

WPStIM 

■ 0.0, 

0.0. 

0.0, 

0.0, 

O.C, 

0.0, 

0.0, 

0.0, 

O.C' 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0' 

0.0 ♦ 

0.0' 

o.o. 


0.0, 

0 . 0 , 

0.0, 

0.0, 

O.C, 

O.C, 

0.0, 

0.0, 

O.C' 

o.o. 


0.0, 

O.C, 

c.o, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0' 

0.0' 


0 . 0 , 

0.0, 

0.0, 

0.0, 

0.0, 

O.C, 

0 . 0 , 

0. 0 , 

O.C' 

o.o. 


o.ie*oi» 

0.0, 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 


tLai i»n ,1 MM-lri'IW. nil • Mi'i n Ill ii ii ll. ...I i 







, MCPl 

a 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0* 
0, 


NCPl 

a 

o» 


T1CPP 

a 

0.0, 0.0, 0.0, 0.0, 0.C, 0.0, 0.0, 0.0, 0.0, o.o. 


CPPTL 

a 

0.3E+00, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, o.o. 


MCP2 

a 

0, 


NCP2 

a 

0, 


T2CPP 

a 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 


CPPT2 

a 

0.3E+00, 0.0, 0.0, 0.0, 0.0, o.o, 0.0, 0.0, 0.0, 0.0* 

%% 

j 

MCP3 

a 

0, 

11 

NCP3 

a 

0, 


T3CPP 

a 

0.0, o.o, 0.0, 0.0, o.o, 0.0, 0.0, 0.0, 0.0, 0.0, 

1 ? 

CPPT3 

a 

0.3E+00, o.o, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

il 

MCP4 

a 

0, 

8a 

NCP4 

a 

0, 


T4CPP 

a 

0.0, 0.0. o.o, 0.0, O.C, 0.0, 0.0, 0.0, O.C, 0.0, 


CPP14 

a 

0. 3E+00, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 


ACP3 

a 

0, 


T5CPP 

a 

0.0, 0.0, 0.0, 0.0, O.C, 0.0, 0.0, 0.0, 0.0, 0.0, 


CPPT5 

a 

0.3E+00, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 


MCP< 

a 

0, 


i> 

M NCP6 

a 

0, 




s> 

T6CPP ■ 0. 0» O.C, 0.0, 0. O f O.Of O.C, OtOf 0.0, 0.0, 0.0, 

CPPT 6 * 0.3C+00, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

NCP5 • 0 , 

SEND • 

PVftCLVSIS CAS KINETICS INPUT 


NUMBCft OP EQUATIONS 20 


EXPi .ENTS OF PATE EQUATIONS 

1211111 112111111001000000 
000000111010011111100000C 

ooooooooooooooooooooooooo 
1111002221000112001000000 
101100 l 1 1 * 0001 10001000000 

00000000000 0 0010001000 000 


# 


lyw w w ^ iw iw^i *! " * mwu.mn mmm w di i mu ' i «t i wwmPMuji 


AF 

SF 

AEF 

i 

7.600000006*14 

0. 

1. 03000000E*02 

i 

l.COOOOOCOE-32 

0. 

0 . 


1.5COOOOCOE*11 
1. OOOOOOOOE- 3 2 

0. 

0. 

1 .62000000E *C 1 
1.00000000E*00 

2 

3. 14000000E*15 
l. OOOOOOOOE— 32 

0. 

0. 

6.010000006*01 

1.000000006*00 

3 

2.570000006*08 

0. 

2. 3000C0C0E*0l 

4 

1.CCC000C0E-3O 

0. 

0. 


2.14C000006*10 

0. 

4.000000006*01 

5 

l.COOCCOCOE-30 

0. 

0. 


1.400000006*1 L 

0. 

5.200000006*01 

6 

l.COCCOOCCE-30 

0. 

0. 


2. 400000006*19 

8.600000006-01 

1 .032000 00E *02 

7 

2.00000000E*18 

1.0000000CE*00 

0. 


2.800000006*21 

1,5000CCCCE*00 

2.250000006*02 

8 

1.600000006*20 

1 . 5O0CCOOCE *00 

0. 


2.500000006*16 

5.OOCOCCO0E-01 

1. 180000006*04 

9 

2.200000006*13 

1.000000006*00 

0. 


3.000000006*12 

-5. OOOOOOOOE -01 

9.92CC0000E*01 

10 

3.000000006*13 

0. 

6.000000006*03 


3.000000006*12 

0. 

4.800000006*00 

11 

1 • OOOOOOOOE- 30 

0. 

0. 


2. 860 000 COE* 06 
1* OOOOOOOOE- 30 

0. 

o* 

6.080000006*01 
1. 80CCC0006*00 

12 

8.780000006 *06 

0. 

7. 71000000E*01 
0. 

13 

1.CCCG00O0E-30 

0. 


7 .8000CGC0E* l 1 

0. 

4.540000006*01 

14 

2.100000006*12 

0. 

0 • 


l.CCC0C000E*20 

0. 

5# OOOOOOOOE *“0 1 

15 


1 . 000000 00E+ 17 

0. 

0. 


2.CCOOOOOOE+14 

0. 

2. 800000 OOE+Ol 

16 

2 .60 0000 00 E* 12 

0. 

6* 3 8 0000 OOE +01 


1.20000000E+12 

-i.ocoooooce+oo 

8.50C00000E+01 

17 

1 . COOOOOOC E-30 

0. 

0. 


9.255000GGE+03 

0. 

7. 00 0000 OOE+Ol 

18 

1 . COOOOOOOE-30 

0. 

0. 


1 • 700000C0E +22 

1. 30000000E+00 

1. 18000000E +02 

19 

4.5OCCC0C0E+21 

1. 5 00 0000 CE +00 

0. 


6.9C0000C0E+13 

0. 

i.eoooooooE+oi 

20 

7.6C0Q0000E+12 

0. 

l.OOOOOOOOE+CO 



NG-23 

SI 


NS- 1 


MM- 4 


XFM 


1.363250C0E+00 
-7. 12450000E-01 
3.04 J69C0CE + 00 
2* 846C8 50CE + 00 
2.8545760CE+00 
3.69161500E+00 
3. 59 761 3CCE +00 
3.7189950CE+00 


1. 85 605 000 E -03 
7.3406500CE-03 
6.U870000E-C4 
4 • 19320000E-C3 
1.5976000CE-03 
-1.33330000E-03 
7.8145600CE-C4 
-2.51670000E-03 


*7. 667 500 00E— 07 
-5. 5262 00 C0E-06 
-7.39900000E-09 
*9» 611 90000E-06 
•6.256tOCOOE-07 
2.65031000E-06 
■2. 238 7C0 C0E-07 
8.5837400dE-06 


1 .5104000 0E-10 
1*51 400000E-09 
-2.03000000E-11 
9.51230000E-09 
1. 13 160 00 OE -10 
-9.76900000E-10 
4.24900000E-11 
-6.29990000E-09 


-1. 13900000E-14 
-2.38200000E-14 
2 .45900000E-15 
-3.30900000E-12 
-7.69000000E-L5 
-9.97700000E-14 
-3.34600000E-15 
2. 7082 0000 E- 12 


l.OOOCE+OO 

H 

0 

5.2098E+04 

2.000CE+00 

M2 

0 

0. 

1.5000E+01 

CH3 

0 

3.3491E+04 

1.600CE+01 

CH4 

0 

•1.7889E+04 

2.600CE+01 

C2H2 

0 

5.4194E+04 

2.8CCC6+01 

C2H4 

0 

1 * 2496E + 04 

3.00001^01 

C2H6 

0 

-2.031 7E+04 

7.80Ci'E + 01 

C6H6 

0 

1.9779E+04 

2 .BQCCE+Ol 

N2 

0 

0. 

1.400CE+01 

N 

0 

1 . 1302E+05 

1.700CE+01 

NH3 

0 

— 1 • 1040E +04 

3.0GCCE+ 31 

NO 

0 

2.1600E+04 

4.6OC0E+31 

N02 

0 

8.0075E+03 

6.2JC0E+01 

NC3 

0 

8.0075E+03 

4.40C0E+31 

N20 

0 

1 .9490C+04 

2.70C0E+01 

HCN 

0 

3.11856 +04 

1.600CE+01 

0 

0 

5.9557E+04 

3 .20 CCE + 31 

02 

0 

0 • 

4.80CCE+01 

C3 

0 

5. 20986^04 


et-i-oooooooooooonj'OjiOHooo 



i.eoooE*oi 

H20 

0 

-5.77986*04 

1. 70006*01 

OH 

0 

9.31256*03 

4.40006*01 

CC2 

0 

~9»4052E*04 

2.8CC06*01 

CO 

0 

-2.64166*04 

1 .20006+01 

c 

1 

0. 


CCHPOSIT ION OF GASES PRODUCED BV NPVR PYROLYSIS REACTIONS 


.0400E-01 

.20006-02 

•OOOOE-03 

.24006-01 


.76006-01 

.98006-01 


- 0 . 

- 0 . 

- 0 . 

1.0400E-01 

- 0 . 

9.20006-02 

6.000C6-03 

4.240CE-01 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

-c. 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

L .76 006-01 
1.96006- 01 
- 0 . 


- 0 . 

- 0 . 

- 0 . 

1.04006-01 

-C. 

9.20006-02 

6.00006-03 

4.2400E-01 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

-C. 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

1.76006-01 

1.9800E-01 

- 0 . 


- 0 . 

- 0 . 

- 0 . 

1.04006-01 

- 0 . 

9.20006-02 
6. 0000E-03 
4. 2400E-01 
- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

1.76006-01 

1.98006-01 

- 0 . 


— 0 . 

- 0 . 

- 0 . 

1.04006-01 

- 0 . 

9.20006-02 

6.00006-03 

4.24006-01 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

1.76006-01 

1.98006-01 

- 0 . 


- 0 . 

- 0 . 

“ 0 . 

1.04006-01 

- 0 . 

9.20006-02 

6.0000E-03 

4.24006-01 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

-0. 

- 0 . 

- 0 . 

1.7600E-01 

1. 98006-01 

- 0 . 



0. 

2.500006*00 

2.500006*00 

0. 

3. 043696*00 
2.846096*00 
1.00000E*00 
2.80277E*C0 
3.395516*00 
1.000006*00 
l • 1 80006 *-00 
4.245776*00 
2.00000E*00 
4.496606*00 
7.903336-01 
2.C0G0CE>G0 
3.502306*00 
1* 12024E *00 
2.000006*00 
6.391546*00 
6.391546*00 
6.00000E*C0 
7.54906E*C0 
7.549066*00 
0. 

2.854586*00 

3.691616*00 

0. 

2.442236*00 

2.514796*00 

C. 

2.14940E*00 
3. 77162E*00 

0. 

3.152546*00 

4.146956*00 

0. 

4.612936*00 

3.434466*00 

C. 

4.626556*00 

2.382126*00 

0. 

4.626556*00 

2.382126*00 


1 .000006*00 

0 . 

0. 

2.00000E *00 
6. 118706-04 
4. 193206-03 
3.000006*00 
6.250456-03 
4.267835-03 
4.000006*00 
1.095006-02 
-6.912706-03 
2 .OOOCCE *00 
5.26980E-03 
2.34660E-02 
4. OOOCCE *00 
1.1 59C0E-O2 
1.39060E-02 
6.000006*00 
- 0 . 

- 0 . 

6.000006*00 

- 0 . 

- 0 . 

0. 

1.597606-03 

-1.3333C6-03 

0. 

1.227616-04 
-1.12430E-04 
3.C00C06 *00 
6.492856-03 
-4.862106-04 
0. 

1.460006-03 

-4.11970E-03 

0. 

2.63866E-03 
2.223 4 36-03 
0. 

3.021686-03 

1.035C5E-02 

0. 

3.021686-03 

1.035056-02 


0. 

0. 

0. 

0. 

-7.399C06-09 

-9.611906-06 

0. 

-2.2 8920E-C6 
2.023276-07 
0 . 

-4. C62C06-06 
3. 160206-05 

0 . 

-1.840006-06 

-3.55420E-05 

0 . 

-4.474506-06 

2.65680E-06 

0 . 

-C. 

- 0 . 

0 . 

- 0 . 

—0 . 

2 . OCOOOE *00 
-6 . 25660E— 07 
2.65031E-06 
1 . C0O0OE*OO 
— 8. 45920E-C8 
2.564756-07 
1.000006*00 
-2.26550E-C6 
9. 87420E-06 
1 • CCCCOE* 00 
-5. 7C70O6-C7 
9.692246-06 
1 . CC0C0E*00 
-1.094906-06 
6. 714906-06 
l.CCCCOE*CO 
-1.215606-06 
-1.116706-05 
2.000006*00 
-1.2 15606-06 
-1.1 16 70E-05 


0. 

0. 

0. 

0. 

-2.03000E-11 

9.512306-09 

0. 

3.899306-10 
-1. 15480E-09 
0. 

7.137006-10 

-2.571506-08 

0. 

3.105406-10 
2. 75510E-08 
0. 

7.94500E-10 

-1.156006-08 

0. 

- 0 . 

- 0 . 

0. 

- 0 . 

- 0 . 

0. 

1. 13160E-10 
-9. 76900E— 10 
0. 

2.140106-11 

-3.24600E-10 

0. 

3. 73940E-10 
-9. 567906-09 
1 . OOOOOE *00 
1. 062806-10 
-7.863406-09 
2 * COOCOE *00 
2.C8130E- 10 
-9.742806-09 
3.000006*00 
2.285606-10 
6.95.M321-C9 
1.000006*00 
2.285606-10 
6.95832 E-09 


0. 

o. 

2.45900E-15 

-3.309006-12 

-2.528006-14 

4.128806-13 

-4.749006-14 

9.510306-12 

-2.00000E-14 

-8.448006-12 

-5.323006-14 

5.238706-12 

- 0 . 

- 0 . 

- 0 . 

- 0 . 

-7.690006-15 

-9.977006-14 

-1.251006-15 

1.259506-13 

-2.361006-14 

3.131306-12 

-7.372006-15 

2.23100E-12 

-1.465006-14 

3.721306-12 

-1.585006-14 

-1.87800E-12 

-1.585006-14 

-1.878006-12 


2.547056*04 

2.547056*04 

-8.549106*02 

-9.672546*02 

1.578756*04 

1.564986*04 

-9. 85560E *03 
-1.018666*04 

2.563706*04 

2.625506*04 

4.543906*03 

5.332896*03 

-1.210176*04 

-1.21017E*04 

7.29743E *03 
7.29743E *03 

-8.901746*02 
-1 .062836*03 

5.614806*04 

5.612786*04 

-6.401966*03 
-6* 7281 OE *03 

9.0522 1E*03 
9.744796*03 

2.340386*03 

2.864776*03 

8.135666*03 
8. 72300E *03 

8.135666*03 

8.723006*03 


-4.600106-01 

-4.600106-01 

-1.648006*00 

-1.411806*00 

5.684126*00 
2. 703756*00 

1.25060E*01 

-9.175506-01 

-3. 144816*00 
1.400506*01 

2.466706*00 

1.583786*01 

- 0 . 

-0. 

- 0 . 

- 0 . 

6.390296*00 

2.287506*00 

4.492576*00 

4.119306*00 

9.238916*00 

1.465406*00 

6.94465E *00 
2.569436*00 

1.367646*00 

8.408466*00 

-1.146366*00 

1.02270E*02 

-1.146366*00 

1.022706*02 




l.COOOOE*OC 

3.65380E*00 

1.000006*00 

3.44363E-03 

l.CC0C0E*00 

-1.25850E-06 

0. 

2.16920E-10 

-1.430006-14 

1.442186*04 

2.37260E*00 

2. 168126*00 

1.07290E-02 

-1. 50880E-05 

1. 19330E— 08 

-3.70000E-12 

1.46829E*C4 

9.281026*00 

0. 

2. 537266*00 

0. 

- 1 . 842206-05 

0. 

-8.8C180E-C9 

1.00000E*00 

5.96430E-12 

-5.57400E-16 

2.923006 *04 

4.946796*00 

3.02189E*00 

-2.173706-03 

3.75422E-06 

-2. 99470E-09 

9.07770E-13 

2.91372E *04 

2.646016*00 

C. 

3.55761E*00 

0. 

7.81456E-04 

0. 

-2 • 238 70E-C7 

2. COCOOE *00 
4.24SOOE-11 

-3. 346006-15 

- 1 • 19279E *03 

o. 

3.719006*00 

-2.5167CE-03 

8.583746-06 

-8.299906-09 

2.70820E-12 

— t .05767 E *03 

0. 

0. 

2.537?6E*CO 

0. 

- 1 . 64220E-05 

0 . 

ecieoE-09 

3.00COOE *00 
5.564306-12 

-5.57400E-16 

2 .92300E *04 

4.946796*00 

3*02 1896*00 

-2.173706-03 

3. 75422E-06 

-2.954706-09 

9. 07770E-13 

2.91372E*04 

2.646016*00 

C. 

2.6TG75E*00 

2.CCQCCE*U0 

3.03170E-03 

C. 

-8. 535C0E-07 

l.OOOOOE *00 
l. 17900E-10 

-6.197006-15 

-2.988906*04 

6.883836*00 

4. 156506*00 

-1. 7244CE-03 

5.69820E-06 

-4. 59300E-09 

1.423406-12 

-3.028886*04 

-6. 861606-01 

A 

2*884556*00 

1 .CQCCCE *00 
9.98350E-04 

0. 

-2. U8C0E-07 

1.000006 *00 
1.4 80 306-1 1 

-3.84500E-16 

3.881186*03 

5.559706*00 

3.823476*00 

-1.118 706-03 

1.24668E-06 

-2. 10300E-10 

-5.25400E-14 

3.58528E *03 

5. 82530E-01 

I .COOCOE *C0 
4.41 290E *00 

0. 

3. 192306-03 

C. 

-1.298C0E-06 

2.000006*00 

2.41500E-10 

-1.674006-14 

-4.094406*04 

-7.28760E-01 

2. 170106 *00 

1 .037806-02 

-1.C7340E-05 

6. 34 592E-09 

-1.62800E-12 

-4. 83526 E *04 

1.066446*01 

1.C3CCC6 *C0 
2.9 51206*00 

0. 

1.5526CC-03 

0 . 

-6. 19110E-07 

1 . OOOOOE *00 
1. 13500E-10 

-7 .788006 15 

-1.423206*04 

6 . 53 140E *00 

3.787136*00 

-2.171006-03 

5.C7573E-06 

-3. 473806-09 

7. 72 170E-13 

-1.436356*04 

2.633556*00 

1 • COCOOE *C0 
1.363256*00 

0. 

1.8560SE-03 

0. 

-7. 66750E-C7 

0. 

1.5 1040E-10 

- 1 . 1 3900E-14 

-6.496706*02 

-7.989006*00 

-7.124406-01 

7.340c 5E-03 

-5. 52620E-06 

1.51400E-09 

-2.38200E-14 

-6.805336*01 

2.793266*00 


(LAST OATA CAKC REAC 


TABULAR LIST CP S6LECTEO OATA EACH KAREi 


CP TABLE SP6CIHC H t AT 
TEMPERATURE 

7.2CCC00006*02 

l.C8COOCCOE*C3 

l.44CCO:;C06*03 

i. acccctcot*cj 


C6 CHAR 
v P 

2. 504000006-01 

2 . •moooooE-ci 

3.4780000CE-01 

3.Sc5C0JCCE-01 




■ - *##**• vm&* 1 *>* <. 




2. 16000000 E *03 
2. 520G0000E *03 
2.574000 00 E* 03 
7.CCCCCQ00E*-03 

CPP TABLE SPECIFIC FEAT 
TEMPERATURE 
0 . 


4. 452000 OOf -01 
4. 939 000 OOE -01 
5. 012000006-01 
5* OOOOOOOOE-Ol 

CF UNCHARRED MATERIAL 
CPP 

3.COOOOOCOE-OI 


CP8AR TAELfc SPECIFIC FEAT OF UASECUS PRODUCTS OF PYROLYSIS 
TEMPERATURE CPBAR 

1.0GC0C0C0E*C3 1 . 00000 CCCE *C 3 

2. 0C0000C0E*03 2. OOOOOOOCE *03 


CM TABLE THERMAL CONDUCTIVITY OF CHAR 


TEMPERATURE 

5.O0CDOOOOE+02 
I * CCCCOOCOE *03 
1.5CCCOOOOF*03 
2.00000000E *0 3 
2* 5C000000 0 03 
3.CCOOOOCOF *03 
3.5CCt:CC0C6»O3 
4* QCCU00C0t*03 


CK 

1.40000000E-05 
l .400000 OOE -05 
3. 68000000E-05 
8.13000000E-C5 
1. 17000COCC-C4 
1.67000000E-C4 
2.C3C00CC0E-C4 
2. 0000 00 COE -04 


CKP TABLE THERMAL CONDUCT 1VITY OF UNChARREO MATERIAL 
TEMPERATURE CKP 

5.000000006*02 I.40000000E-C5 
1.0CC()O0COE*03 1 • 40000000 E-0 5 
1.500000006*03 3. B8000000E-05 


Hh TABLE LOCAL ENTHALPY CF FLUIC AT FALL 


TEMPERATURE 
5. 000000006*02 
1 .000000006*03 
1.50C0C0C06 *03 
2.000000006*03 
2.5CC000C0E*C’ 

3. OCC 000006*03 
3. SCCCOO 006*03 
U 4.000000006*03 

4. 5 C COCO 00 c* 03 
5. 00 COO 0006 *03 
5.5CCCCCC0C*03 
6* CCCt/UOCOL *03 
6.500000006*03 
7.C0C) OOOOE*Oi 
T8 TABLE TEMPERATURE 
TIME 


HH 

-1.000000006*01 1 
1. 11000000E *C2 
2. 390 JOG COE *02 
3.75000000E*02 
5. 17C03000E*02 
6.650000006*02 
8. 100000006*02 
9.800000006*02 
1.1 73000 00 E *03 
1.384000006*03 
1 . 700C COCCE *03 
2. CbO 000006*03 
2.550000006*03 
3.000000006*03 

TO WHICH BACK SURFACE RACUTES 
TB 


0 


C 


hE TABLE LOCAL ENTHALPY EXTERNAL TC BOUNDARY LAYER 
TIME HE 

0. l.92Z000CCE*C4 

HC TABLE HEAT CF SUBLIMATION OF CHAP 
TIME HC 

2. 0CCC000OE*03 1. 20000C OOE *04 



Q{. TABLE CONVECTIVE 
TIME 

0. 

1 *00000000 E *01 
2.occoooooe*oi 

. 3 . 000000 COL ♦01 
4*OCCOGOOOt*01 
5.0CCOG0GGE*01 
6.0CCo00G0t*0l 
7.CCC000C0I *01 
8.ooogoogoe*oi 
•j.ccccoccot *oi 
l.OGOCOOCOL*02 
L * l CCCCO CO fc *-02 
1.2CCC00COE*02 
l.30000000E*02 
1.4CCOOOCOE*02 
1.5CCC0000E*02 
l, •'C OCCOOOl *02 
1 700^20001*02 

1. accooooot *02 

1.9CCCOOOOt*02 

2. CC0000C0E *-02 
2.lCCCG000t*02 
2. 2 00000 00 E>02 
2.3CCOOOOOL*02 
2.4CC0C0C01 *-02 
2. 500 OjOCOL*O 2 
2.6CCCCC CO t *02 
2* 7C0 OOOOOL *U2 
2.8CC00G00fc*02 
2. 9 0000 000 t ♦02 
3.000000000*02 
3* l COCOOCOE*02 
3.2CCCOOCOt*02 

3. J3CCCGC0E*92 
3.4COOOOCOfc *02 
3 . 5CCCCOOOE * 02 
3.6CCOOOOGt*02 
3.7C0C3QCOE*O2 

3. BCuCoOCO E *02 
3.9COCCOOCL ♦ 02 

4. ccrco0(-0f *02 


HEATING RATE 

OC 

3. 65CCQCCCE *C0 
7. 500000 COE ♦ CO 
1.680000COE*01 

4. C60C00CCE *C 1 
9. 61000C OOE *01 
1 . 92 7 COC COE *02 
3. 1 7709000E *02 
4. 3B500000E *02 
4.90100000E*C2 

4. 344 000 OOE *02 
3.3J7000COE*02 
2.61 OOGOCCE *02 
1. 698C000CE*C2 
1. 397000CCE *C2 
1. 16300000E *02 
S.790000COE*C1 
8.31 0000 OCE *01 
7. 120000 00E ♦Cl 
6.1 5( OOt OCt *C 1 

5. 33000000E *01 
4.72COOOCOE*01 
4.22000000E *C 1 
3. (120000 OOE *0 1 
3 .4900 00 OOE *01 
3* 2 1000000E *01 
2.95OC0CCCE*0l 
2. 71000 OCOE *01 
2. 46 00 00 COE *01 
2. 1BOOOOOOE *01 
1.87000000E *01 
1.5400 00 00E*C1 
1. 22CUOCOOE *01 
5.2000000CE*CO 
6.7000COCCE *00 
4. 700000 CCE *00 
3. 300000COFKC 
2«2000000CE*CO 
2.200CCUCCE »C0 
1.4DOCOOOCE*CO 
1. C0000C0CE*00 
6.000000C0E-C 1 




rnatmuummikitiumiUit 





*.iCOOOOOOEM>2 

«.2CC000C0E*C2 

4.3CCCOOCOE*-02 

4. <*CCCOOOO£t02 
A.5CCCOOOOE*02 
*.6CCOCOOOt*C2 
4.7CCCCCC0E »02 
4.8CCCCGOO[>02 
4 • 900000001 ♦•02 
5.O0C000O0E»02 
5. lOOOOOOOt ^02 
5.20000000E>02 
5.3CCOOOOOE»02 
5»<r0C0J000E»02 
5»50000CCOt*02 
5.6OC00000fc*C2 
5.7CCC00C0E+O2 

5. UCCC0000EO2 
5.9CCO0000E^02 

6. C0CU00G0E^02 
6. lCCGOOGOt + C2 
1. 30C00000E»03 
2.000000COE»03 


A.OOOOOOOOE-Ol 

3.00000000E-C1 

2.00000JCOE-01 

1.5COOCOOOE-01 

1. tOOOOOCOE-Cl 
8. OOOJOOOOE-02 
6.00000000E-C2 
5. 033COO0CE-02 
A. 00000000E-Q2 
3.00000000E-02 
3.00000000E-C2 
3. OOCOOOCOE-C2 
2. OOOOCOOOE-02 

2. OOOOOOCOE-02 
2.030000COE-02 
2. 000000 00E-02 
2. OOOOOOCOE-02 
2* OOOOOuOOE-02 
l.tOOOOOOOE-02 
l.OOOOCGOOE-02 
1»OOOOOOOOE-03 
1. 000000 OOE-03 
l„ OOOOOOOOE-03 


QR table raciant seating rate 
TIME QR 


0 . 


5. COOOOOOOE-02 


P-RATIO TABLE 

TIME PRAT 

0 . I .OOOOOOOOE *00 

Q-RAT10 TABLE 

TIME QRAT 

0. 1.0000000CE*00 


PRESSURE TABLE 

ume PRESSURE 

0. 3.88300000E-C3 

5.COCOOOOOE»00 7.9A500000E-03 


ALPHA TABLE 
TIME 

0. 


ABSORPTIVITY OF CHAR SURFACE 
ALPHA 

1.19900000E*00 


ALPHA1 TABLE 
TIME 

0. 


ALPHA1 

5.00000000E-01 


ALPHA2 TABLE 
TIME 


ALPHA2 



1.000000006*00 


MCDOT TABLE RATE 06 CHAR LOSS-OPTICA 2 
TIME MCDOT 

0 . 0 . 


X TABLE CHAR THICKNESS- 
TIME 

0. 

DELTAU TABLE 
TIME 

0. 

b.ocooooooe*oo 

1 . 00 CO 00 006*01 
4 . 000000006*02 

HCCMB table 


TEMPERATURE 

l.8CCCOOOOE*03 
2. TCCC00C06»03 
3.600000006*0 A 
4. 9CCC 00006*03 
5.403000006*03 
6.30000CU>t *03 
?.2CC0G3C0t*G3 

SOP ■1.6350306-01 
Beta • o.ocoo 
TC • 0. 

EIPJ * 0. 

AEXP - 1.000000006*10 

QR06R- 1.0 
ABEXP- 0. 

I • 26 

L DISTANCES 

2 . 0 CCCCCC 06-02 
2 .CC 0000006-02 
2 . 0 CGCCOC 06-02 
2 . COO 000006-02 
2. 000 000 00 E— 02 
2.000000 C0E-02 
2 . 0 C 000 C 0 O 6-02 
2 .OCcCCv'Cwt — 02 
2 . C C C 3 c o v w i — C 2 
6 .CCo 3 dw ; _-Ll 
4. C C60u0c -cl 
6 . CC 60 CCE- C l 


CPTION 3 
X 
C. 


OELTAU 

T.B 125000 CE-C 3 
1 . 5625 C 00 CE-C 2 
6 . 25000 COOE -02 
1 . 000000006*00 

PRESSURE 

1.000 000 00E-01 


4 . 100000006*03 
4 . 264000006*03 
4 . 454000006*03 
4 . 87100090 E *03 
6 . 26 5 COO JOE *03 
1 . C 2200 ) 00 E *04 
1 . 354 O 0 COO 6 »C 4 


1.000000006*00 


4 . 110000006*03 
4 . 266000006*03 
4 . 4 4400 0 00 E * 03 
4 . 69 TOOOOOE *03 
5 . 295 CCC 006 *0 3 
6 . 995000006*03 
1 . 305000006*04 


1.000000006*01 

4 . 1 1400000 E *03 
4 . 26600000 E » 03 
4 . 446000006*03 
4 . 656000006*03 
4. 904 00 OOOE* 03 
5 . 679000006*03 
7 . 1 34000 OOE * 03 


PO • 0.0000 ROP - 
ETA ■ .6000 REFF - 
DEIMP* 0 . 06 LHF- 
EUN » 0. C IP J • 
OEXP • 7 . 650000006*04 APEXP* 
CE - 2 . 32000000 E -01 ELAM . 
BBEXP* 0 . IOPT • 
J • 23 M 


0.0000 

1.00000000 


1 . 240000006*05 

0. 

1 

0 


ROOP • 
COEP ■ 
EPSONE* 
CIJM > 
BPEXP ■ 
HCOMB • 
I 0 PT 1 * 
ERROR* 


1.000000006*02 


4 . 110000006*03 
4 . 266000006*03 
4 . 445000006*03 
4 . 643000006*03 
4 . 884000035 *03 
5 . 245000006*03 
5 . 669000006*03 

5.5000 

2 . 024000006*01 

7 . 500000006-01 

0. 

1 . 960 000 OOE *04 

0. 

1 

1 . 000000006 -03 


2 . 30000300 E -32 
2 . 000000006-02 
2.00000000E-02 
2 . 000000 C 06-02 
2 . 000000 OOE -02 
2.00000000E-02 
2 . 000 d 30 i'OE-C 2 
2.0:0)30':OE-02 
2 .C 0 O 0 JJ 0 CE-C 2 

.0543 

. 854 C J JE -01 

. 854000 E -01 


2 . 0 COOOCOOE -02 

2.ooooooccr-o2 

2 . 000000 COE - 02 
2 . 000000 OOE - 02 
2. 000000 OOE- 02 
2 . 000000 C 06-02 
2 . 0000 OC COE- C 2 
2 . 0 CO 30 CC 06-02 


2 . 000000006-02 

2 . 000000006-02 

2.00000000E-02 

2 . 000000006-02 

2 . 000000006-02 

2.00000000E-02 

2.00000000E-02 

2 . 000000006-02 


2.00000000E-02 
2 . 000000006-02 
2 . 000000006-02 
2 . 000000006-02 
2 . 000000006-02 
2 . 000000006- 02 
2 . 000000006-02 
2 . 000000036-02 


2 . 000000006-02 
2.00000000E-02 
2 . 00030000 E -02 
2 . 000000006-02 
2.00000000 E-02 
2 . 003000006-02 
2 . 000000006-02 
2.00000000E-02 


1.0164006 » C l 
1 . 0164006*01 
1 . 0 1640 C 6 * C 1 


0 . 0025006-01 
B. 0025006-01 
B.C 02500 E -01 


1 . 3612506*00 
1.36 12506 * 00 
1 . 36 I 250 E *00 


6 . 0805006*00 
6 . 0885006*00 
6 . 0885006*00 


1 . 3203006*01 

1 . 3200006*01 

1 . 3209006*01 


n o w. ‘ 






«.6UOOOE-Ot 
8.CC6000E-01 
A.CS6000E-01 
8.0060O0E-01 

6.CC6G00E-C l 

6.CC6 JOOE-Ul 

6.CC6000E-01 

6.CC6GGCE-01 
6.0Ct000E-Cl 
fc.CCGOOGE-Cl 
6 . UC4 COOE-C 1 
6.CC60u0E-0| 
6. OCtCGGE-Ol 
6. COCO 00 =- 01 
6.CC6CCCI-01 
6.CC60C0E-01 
6.CC6uOOE-OI 
e.ctbOudE-ci 
4.C060G0E-C1 
6.CCSOJ0E-01 
6.C060006-GI 
6. CC6QG06-C l 
6. CC 6 O 0 OE-CI 
6.CC6G006-CI 
6. CC61 OOE-Cl 
*.CCg0006-OI 
6.CC60O0E-C1 
6 .OO 6 GOOE-OI 
6.CC6000E-G1 
6 . CG60U0E-C 1 
6 . C C630Q t-0 l 
6.CC60ugE-C1 
4.0C6000E-G1 
6. CC600GE-G 1 
4. QC600C6-01 
4.CC6CC06-C1 
6. CC6GGCE -0 1 
6. JC6C0CE-GI 
6.CCtOOQt-Gl 
6 .CC 6 O 0 OE-CI 
6. CC6U00E-C l 
6.0060006-01 
6. 0C60006-01 
«. CC60G06-01 
0.CC6O<.'Jt-Cl 
AcCCcwOG'-'Cl 
b.rc^w-'GL-tl 
i.OCCCJOE-Ol 
O.CC600 Jt-01 
6.CC6GOOE-01 
6.CC6G006-01 


7.8540006-01 
7.654000E-01 

7. 8540006-01 
7.6540006-01 

7.8540006-01 

7.8540006-01 

7 . 0540006- 01 
7.6S40OJE-01 
7.054COOE- 01 

7.8540006- 01 
7.6540006-01 

7.0 540006- 01 
7.854COOE- 01 

7 .8540006 - 01 

7.8540006- 01 

7.0540006- 01 

7.e54O00E- 01 
?.654000t-0l 
7.854000E- 01 
7.85400OE- 01 
7.85 40 00 £- 01 
7.8540OJ6- 01 

7.0540006- 01 

7.B54GO0E-01 

7.8540006- 01 

7.U54000E-01 
7.054uOUE- 01 
7 .83 40006- 01 
7.854000E-01 

7. 8540006- 01 
7 . 854000C -0 l 
7 • U5 4ooOE -0 l 
7.854000E-01 
7. 85400GE* 0 1 

7.8540006-01 
7.85400OE-01 
7.854003E-01 
7 • 0540006- 01 
?.854uOoE-Ol 
7.854000E- 01 

7.8540006-01 
7.854000E-01 
7.854000E-01 
7. 8S4000E - 01 
7.854000E-01 
7. , J54„JU'‘- 01 
7.e54« oot- 01 

7. 0J.OOOE- 01 
7.6 54C 306* 0 1 

7. 8540006-01 
7.854OOOE-01 


1*0164006401 
1.0164006*01 
1*016 400 E 4 C 1 
1 . 01 6400E *01 
l. 0 1 640 u£ ♦ C l 
1*016400 E*C1 
1* 01640GE*01 
1 . 01 640C E ♦ 01 
1.016 40 06*01 
1*01 6400 E * C 1 
1*01 64 00E 4 0 1 
1*0164 00 E401 
1.0164006*01 
1* 0164006 *01 
1 ,0164306*0 1 
1*01 64GOE *01 
1 • C 164GC6 *01 
1.0164006*01 
1.0 1 640GE 4 0 1 
1.0164006*01 
l.C16400'£»01 
1.0 1 6 40ub * 0 1 
l.C16400E»01 
1. C 1 6400E 40 1 
l. 0164006*01 
1*0 164CGE 40 1 
1.0164006*01 
l* C 16400E 4 01 
1,0 16400E 4 0 1 
1.0164006*01 
1. C 164G0E 401 
1.0164006*01 
1*01 6400E >01 
l.016400E»01 
1* 01 6400E 40 1 
1*0164006401 
I.C1640CE401 
1.0164006*01 
1*01 64Q0E* 01 
1.016400E*0l 
1 . 01 640oE *0 1 
1. 0164006*01 
1.C16400E401 
1*01 64 00E * 01 
1.0164006*01 
1.01 640oE*01 
l.Ct64uOE*01 
1. 0164006*01 
1 . 01 6400E *01 
I.C16400E*01 
l. 0164006*01 


8. 0025006-01 
8* 0025006-01 

8. C0.1500E-01 
8.oo;:5ooe-oi 

8. CC25006-01 
6. 002S00E-01 

8.C02500E-01 

8.0025006-01 

8.0025006-01 

8.C025006-01 

8.0025006-01 
8* C02 5006-01 
8 • 002 5006-01 
8. C02500E-01 
B.002500E-81 
8* 0025006-01 
8.CC2500E-01 
8.0Q2500E-01 
8. C 02 5006- 0 1 
8.C02500E-01 
8. C32500E-01 
8. 002500E-01 
0.OO25OOE-O1 
8.C02500E-01 
6* 002 50GE-01 
a. 002 5 00E- 0 1 
8. C02 500E-01 
0. 002 5 OOE- 01 
8.C02000E-01 
8.002500F-01 
8. C 02 5 OOE-Ol 
0.OJ25OOE-OI 
8.C02500E-0I 
8 • C02 5CCE-01 
8.002500E-01 
8.002500E-01 
8.002500E-01 
8.C02500E-01 
8 • 002500E- 01 
8* 002 SOOE-Ol 
8.002500E-01 
8. 002 500E-01 
B.002 r C06-01 
8. 002 500 E-Ol 
8.C02500E-01 
8 « 002 500':— 0 1 
8. C02500E-01 
8* CO? 5006-01 
8* C025006-01 
8.0O2500E-01 
8. 002500E-01 


1.3*12506*00 
t. 361250E *00 
1* 36 12 50E* 00 
1.361250E400 
1.3612506*00 
1.3612506*00 
1. 36 1250E*00 
1. 3612506*00 
1.361250E*00 
1. 36 1250E » 00 
1-3612506*00 
1. 36 1250E * 00 
1. 36 1250E * 00 
l. 3612506*00 
1* 36 12 50E * 00 
1 . 36 1250E *00 
1. 3612506*00 
1.3612506*00 
1.3612506*00 
1. 36 1 250E ♦ 00 
1. 3612506*00 
1.3612506*00 
1.3612506*00 
1.3612506*00 
1.3612506*00 
1.36 1253E*00 
1. 36 1 2 50 E * 00 
t. 3612506*00 
1.3612506*00 
1.3612506*00 
1. 36 1 2 50F * 00 
1.361250E*00 
1.3612506*00 
I. 36 1250E * 00 
1.3612506*00 
1. 36 12 50E *00 
1.3612506*00 
1.36 1 253E * 00 
1. 36 1250E * 00 
l. 3612506*00 
1.36 1^506 4 00 
1.3612506*00 
1.3612506*00 
1.3612506*00 
1.3612506*00 
1. 35 12 50 f *00 
1.3612506*00 
1.3612506*00 
1.3612506*00 
1.3612506*00 
1. 36 1250E *00 


6.0885006400 

6.0885006400 
6 . 388500E »00 
6 . 088500E *00 

6.0885006400 
6. 0885006*00 
6.0885006*00 
6. 0885006*00 
6. 0805006*00 
6. C88300E*00 
6. 0885006 *00 
6. 088503 E *00 
6.0885006*00 
6. 0885006*00 
6. 0885006*00 
6. 0885006*00 
6. 0883006*00 
6. 088500E *00 
6 • 080 5 0 OE *00 
6. 0685006*00 
6* 08850 JE *00 
6. 0885006*30 
6.0685006*00 
6. 0883006*00 
6. 088500E *00 
6.0885006*00 
b. 088S00E *00 
6.0885006*00 
6. 0885006*00 
6. 0865005*00 
6.0383036*00 
6 . O083OOF *00 
6 . 0385006 *00 
6. 0B8500E *00 
6.0885006*00 
6.0805006*00 
6. 0805006 *00 
6. 088509E *00 
6. 088500E *00 
6.0805006*00 
6. 0885006*00 
6.0885006*00 
6. 0885006*00 
6. 0885005*00 
6. 3885006*00 
6. 0685006 *00 
6.086333E*00 
6.O803OOE*OO 
6. 0035006 *00 
6. 0065006*00 
6. 0885006*00 


• 3200006*01 
.3200006*01 
.3200006*01 
.3200006*01 
.3200006*01 
.3200006*01 
. 320000E*01 
.3200006*01 
.3200006*01 
.3200006*01 
. 3200006*01 
.3203006*01 
. 3200006*01 
. 3200036*01 
.3203006*01 
. 3200006*01 

• 3200006 *01 
. 3200006*01 
.3200006*01 
.3200006*01 
. 3203306*01 
.3233006*01 
. 3200006*01 
. 3203006*01 
. 320000E *01 
.3203006*01 
. 3200006*01 
. 3200006*01 
. 320000E 401 
. 3200006*01 
.3200006 *01 
. 320000E*m 
. 32 00 OOE *0 i 
.3200006*01 
. 3203 00E *01 
. 3200006*01 
. 3200006*01 
.3200006*01 
.3200006*01 
. 3200006*01 

• 3200006*01 
. 3200006*01 

• 320000F *01 
l. 320000E *01 
l. 320300E*01 
.. 320100601 
l. 3200006*01 
1.3200006*01 
1.3200006*01 
1. 3200006*01 
l. 3200006*01 



w *m 


m s wum « f % w . »r > ??! >' ivtwiwi i 


4.QC4000E-CI 

4.0C6000E-01 

4.CC6C00E-01 

6.0C6C006-01 

O.CCfcOOOE-Cl 

A.CCb000c-01 

6.CC6000E-U1 


7. *340006-01 
7*6540006-01 
7. 8540006-01 
7.8540006-01 
7.6540006-01 
7.6540006-01 
7.654000E-Q1 


,1.0144006401 
1.0164006*01 
I.C1640QE*C1 
1.0164006*01 
l. 0164006*01 
, 1.0164006*01 
' 1.0164006*01 


а. 0025006-41 

0. 002500 E-Ol 
8. CC2500E-01 

б. 0025006-01 
8.0025006-01 
8. C02 5006-61 
8.0325006-01 


1.3412906*00 
1. 3412506*00 
1.3412306*00 
1. 3612506*00 
1.3612506*00 
1.3612506*00 
1.3612506*00 


4.0683006*00 
6.0885006*00 
6.0885006*00 
6.0885006*00 
6. 0885006*00 
6.0885006*09 
6.0885006*00 


1.3200 006*41 
1.3200006*01 
1. 3200005*01 
1.3200006*01 
1.3200006*01 
1.3200006*01 
1. 320000E*01 






t. 


I 


"V atjw : « tmcr **&**«**• 


j>+ v4u»>:. ; jpv***»*m. x**.**tz- • n w. \nx+*W } 





,007813 






-VTM '•**»:»■ 


►, ym*jr 'K -' -; . i*r*i «* .y iy. «■»*•% • 


OELTAU* o«iC?0U50 

MMPOOT * K 4 C 

0 , •OflOO , 

0 , .COCO / 

0 . .0000 

o, »0000 

0 . *0000 

o. • ocoo 

o. ,coco 

0 , ,0000 

o. .coco 

0 . ,0000 

Q # « cooo 

0 . ,0000 

o. ,cooo 

0 # ,occo 

o. ,0000 

o, ,0000 

o. .0000 

o. ,0000 

o. .ocoo 

o. .0000 

o. .occo 

o. .oooo 

o. . ccoo 

o. .cooo 

o. .0000 

o. .ccoo 

Q # .0000 

o. ,0000 

q m .cooo 

Q a , 0 c oc 

o. .0000 

o. .0000 

o. .ccoo 

o. .0000 

Q, .ocoo 

o. *0000 

o. .cooo 

0 , .0000 

o. .ccoo 


RHO 


.6006 

.7094 

.6006 

.7094 

.6006 

.7054 

.6006 

.7094 

.6006 

.7054 

, 6 G 06 

.7054 

.6006 

,7054 

, 6 G 06 

.7654 

,6006 

.7054 

.6006 

.7854 

.6006 

.7054 

.6006 

.7854 

.6006 

,7054 

.6006 

,7854 

.6006 

.7854 

.6006 

.7854 

.6006 

.7854 

.6006 

.7054 

.6006 

.7054 

,.6006 

.70 54 

.6006 

.7054 

.6006 

.7654 

.6006 

.7054 

.6006 

.7054 

.6006 

.7854 

.6006 

.7054 

.6006 

.7054 

• 6 C 06 

.7854 

.6006 

.7854 

.6006 

.7854 

• 6006 

♦ 7654 

.6006 

.7854 

• 6006 

.7854 

.6006 

.7854 

.6006 

.7854 

• 6006 

,7854 

.6006 

.7854 

.6006 

. 7 e 54 

.6006 

.7854 


10 . 16*0 

10. 1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10. 1640 

10.1640 

10.1640 

10.1640 

10.1640 

10. 1640 

10. 1640 

10. 1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10. 1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 

10.1640 


•8002 

1.9612 

.8002 

1.3612 

,0002 

1.3612 

.0002 

1.3612 

*8002 

1.3612 

.800 2 

1.3612 

.6002 

1.3612 

.8002 

1.3612 

.8002 

1.3612 

• 6002 

1.3612 

.8002 

1.3612 

.8002 

1.3612 

.800 2 

1.3612 

.8002 

1.3612 

• 8002 

1.3612 

.800 2 

1.3612 

.8002 

1.3612 

.0002 

..3612 

.0002 

1.3612 

,800 2 

1.3612 

.8002 

1.3612 

.8002 

1.3612 

.8002 

1.3612 

.0002 

1.3612 

.8002 

1.3612 

.0002 

1.3612 

.8002 

1.3612 

,8002 

1.3612 

,8002 

1.3612 

.8002 

1.3612 

,8002 

1.3612 

.8002 

1.3612 

.8002 

1.3612 

• 8002 

1.3612 

.0002 

1.3612 

.800 2 

1.3612 

• 8002 

1.3612 

• 8002 

1.3612 

.8002 

1.3612 


6.0009 

13.2000 

6.0009 

13.2000 

6.0809 

13.2000 

6.0889 

13.2000 

6.0889 

13.2000 

6.0989 

13.2000 

6.0889 

13.2000 

6.0899 

13.2000 

6.0989 

13.2000 

6.0889 

13.2000 

6.0895 

13.2000 

6.0889 

13.2000 

6.0889 

13.2000 

6.0885 

13.2000 

6.0885 

13.2000 

6.0385 

13.2000 

6.0005 

13.2000 

6. 0905 

13.2000 

6.0005 

13.2000 

6.0085 

13.2000 

6.0885 

13.2003 

6.0085 

13.2000 

6.0095 

13 . 2 C 00 

6.0885 

13.2000 

6.0889 

13.2000 

A . 0(89 

11.1000 

6.0089 

13.2000 

6.0885 

13.7000 

6.0889 

13.2000 

6.0885 

13.2000 

6.0989 

13.2000 

6.0885 

13. 2000 

6.0085 

13.2000 

6.0385 

13.2000 

6.0885 

13.2000 

6.0885 

13.2000 

6.0885 

13.2000 

6.0805 

13.2000 

6.0085 

13.2000 










5J 


UNI 

TUN! 


540.00 

0. 

.0000 

• 600b 

.7854 

10.1640 

• 8002 

1.3612,' 

6.0885 

13.2000 

540.00 

0. 

.0000 

.6006 

.7854 

10.1640 

.8002 

1.3612 

6.0885 

13.2000 

540.00 

0. 

.0000 

.6006 

.7854 

10.1640 

.8002 

1.3612- 

6.0885 

13.2000 

540.00 

0. 
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TABLE I 


DEFINITION OF SYMBOLS 


The following table contains definitions of important symbols used 
in the program. The first section is in alphabetical order by physi- 
cal symbol. The second section contains those items requiring 
definition but not assigned a physical symbol. The latter section 
is in alphabetical order by FORTRAN symbol. 


SECTION I 


Physical 


FORTRAN 


Symbol 

Number 

Symbol 

Definition 

A 

10 

AEXP 

Pre-exponential factor of Arrhenius 
expression for calculating MCD0T. 
Corresponds to specific reaction 
rate. 

A’ 

20 

AP 

Pre-exponential factors of Arrhenius 
expression for calculating pyrolysis 
rate of each pyrolysis reaction. 

B 

40 

BEXP 

Power of the exponential term in 
the Arrhenius expression for MCD0T. 
Corresponds to activation energy 
divided by the gas constant. Calcu- 
lates K in equation (15) as follows: 

K « AEXPe" BEXP/T(1) 

Also used for oxidation calculation. 

B' 

50 

BP 

Pre-exponential factors of Arrhenius 
expression for calculating pyrolysis 
rate of each pyrolysis reaction. 

C 

80 

— 

Concentrations of a specie in the 
pyrolysis gas. 

C e 

90 

CE 

Oxygen concentration in the external 
flow. 

C P 

91 

CPTAB 

Specific heat of char as a function 
of temperature. 


62 


Physical 


FORTRAN 


Symbol 

Number 

Symbol 

Definition 

C ' 

P 

N 

100 

CPPT1 

CPPT2 

CPPT7 

Specific heat of each specie in abla- 
tion material as a function of temper- 
ature 

C " 
P 

no 

CPDPTB 

Specific heat of insulation layer 
as a function of temperature 

C 

P 

120 

CP BAR 

Enthalpy of pyrolysis gas as a 
function of temperature. 

C i+j 

130 

CIPJ 

Concentrated heat sink at back of abla- 
tion material 

C i+j+m 

140 

CIJH 

Concentrated heat sink at back of 
insulation layer 

AE 

150 

— 

Activation energy. Appears in computer 
program only in the form AE/R 

F 

160 

FI 

Constant heat generation function per 
unit volume 

G 

180 

(PC0EFF) 2 

The constant used in the hypersonic 
pressure term. See equation (9) in 
reference 1. 

H 

c 

190 

HCTAB 

Latent heat of process which removes 
material from the front surface as 
a function of temperature and 
pressure 

H 

e 

200 

HETAB 

Enthalpy of the gas stream external 
to the boundary layer. Zero point 
in table should correspond to 
temperature on which cold wall 
heating rate is based. 

H 

w 

210 

HWTAB 

Enthalpy of gas stream at the surface 
temperature. Zero point in table should 
correspond to temperature on which 
cold wall heating rate is based. 

*Ah 

P 

220 

DDELHP 

Heat of pyrolysis for each reacting 
specie. 
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Physical 

Syinbol 

Number 

FORTRAN 

Symbol 

Definition 

Ah 

c 

230 

HC0MBTB 

Heat of combustion of char as a 
function of temperature and pressure. 

Ah f 

240 

DELHF 

Heat absorbed by coolant behind third 
layer. 

i 

241 

I 

Number of stations in first layer 
(1 > 26 is recommended) 

j 

242 

J 

Number of stations in second layer 
(J > 25 is recommended) 

k 

250 

CKTAB 

Thermal conductivity of char as 
a function of temperature 

k* 

260 

CKNPTB 

Thermal conductivity of virgin 
ablation material as a function 
of temperature 

k" 

270 

CKNDPT 

Thermal conductivity of insulation 
material as a function of temperature 
and location 

1 

280 

AL 

Non-dimensicnal distance between 
stations 
i+j-1 

£ AL ■ 1.0. If AL(1) ■ 0. then 

n-1 n 

AL(n) are all equal to AL(n) ■ ^ 

M 

290 

XFW 

Molecular weight of each reacting 
specie in pyrolysis gas 

m 

300 

M 

Number of stations in the insulation 
layer (M 5^ 10 is recommended) 

• 

m 

c 

310 

WCCD0T 

Surface mass loss rate by combustion 

m 

P 

320 

WPD0T 

Local mass flow rate of pyrolysis 
products 

n 

330 

X0RDER 

Order of oxidation reaction. Use 


value of 0.5 or 1.0 


Physical 

Symbol 

Number 

FORTRAN 

Symbol 

N 

340 

NG 

XN°* 6 

Le 

350 

ELAKTB 

P 

w 

360 

(EXPRESS) 

q 

aero 

380 

QAER0 

q c 

390 

QCTAB 

q r 

400 

QRTAB 

q s 

410 

QSTAB 

R 

420 

RR 

R 


430 

RATE 

T 

460 

T 

t b 

470 

TBTAB 


Definition 

Number of chemically reacting species 
in pyrolysis gas 

Product of Levis number to the 0.6 
power and X where X is the weight 
of char removed from the outer surface 
per unit weight of oxygen diffusing to 
the surface 

Pressure at the wall as a function of 
time 

Net aerodynamic heating rate to the 
surface after correction for hot wall 
and blocking 

Reference convective heating rate. 
Stagnation values usually used as 
reference. (QRATT specifies ratio 
of local to reference heating rate.) 

Reference radiant heating rate. 
Ahsorptance, ALPHAT, is altered to 
account for ratio of local to reference 
heating rate. 

Stagnation heating rate 

Gas constant or 

Maximum radius of ablating 

body. Required when I0PT-5. 

Effective reaction rate for gas or 
solid species in pyrolysis gas 
reaction 

The initial temperatures of the fin- 
ite difference stations. Stations 
are numbered from the front surface. 

Temperature as a function of time to 
which back surface is radiating 
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Physical 

Symbol 

Number 

FORTRAN 

Symbol 

Definition 

T 

i+j 

480 

TIPJ 

Back surface temperature at which 
coolant becomes active (TNST0P»O.) 
When TNST0P -1, program stops when 
back surface temperature reaches 
TNBAR. 

T i+j+m 



Same as T... when there is 
another ^ material behind 

ablator 


490 

TAU0 

Starting time 

C f 

500 

ENDTAU 

Stop time if TNST0P - 0. If TNST0P is 
1, set ENDTAU very large and pro- 
gram stops when back surface temp- 
erature reaches TNBAR 

U 

510 

— 

Velocity of pyrolysis gas 

V 

520 

VI 

Surface recession velocity 

AW f 

X 

530 

540 

DELWF 

X 

Coolant consumption rate “ back 
surface 

Thickness of the char layer 

X* 

550 

XP 

Thickness of the virgin ablation 
material 

X*' 

560 

XDP 

Thickness of the insulation layer 

y 

570 

— 

Coordinate normal to the surface 

z 

580 

Z 

Location of fixed stations measured 
from the back surface 

a 

590 

ALPHA 

Name assigned internally to absorp- 
tance of the surface. Value of 
ALPHA is obtained by program from 
table entitled ALPHAT 

a 

c 

600 

AL1 

Name assigned internally to factor 
used to correct blocking effectiveness 
of mass loss from outer surface for 
molecular weight effects and/or 


turbulent flow. Values of AL1 
obtained by program from table 
entitled AL1TAB. 
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Physical 

Symbol 

Number 

FORTRAN 

Symbol 

Definition 

a 

P 

610 

AL2 

Name assigned internally to factor used 
to correct blocking effectiveness of 
pyrolysis gas due to molecular weight 
effects and/or turbulent flow. 
Calculated automatically when py- 
rolysis gas reactions kinetics are 
treated in detail (ICPBAR * 2) 

3 

620 

BETA 

Determine if ablation or transpiration 
theory will be used. If ablation theory 
BETA ■ 1; if transpiration theory, 

BETA * 0. In the former case, a 
transpiration factor ETA must be 
specified 

Y l* Y 2 

621 

— 

Approximate geometric coefficients 

£ 1 

630 

EPS0NE 

Emittance of the heated surface 

e i+j 

640 

EIPJ 

Emittance of back surface when no 
insulation layer is present 

£ i+j+m 

650 

EIJM 

Emittance of back surface when an 
insulation layer is present 

n 

660 

ETA 

Transpiration factor. A value of 0.6 
is frequently used for laminar 
flow. A value of .2 is used for 
turbulent flow. 

y 

— 

— 

Viscosity of gas flowing through char 

p 

665 

RH0 

Density of solids that undergo 
pyrolysis at each finite difference 
station 

p c 

670 

RH0C 

Density of char at each finite 
difference station 

p" 

680 

R0DP 

Density of insulation 

At 

690 

DELTAU 

Time step obtained from table entitled 
TAUTAB by zero order interpolation 


700 

PHIETB 

Parameter equal to M^k 


2RTy 
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Physical FORTRAN 

Symbol Number Symbol 


Definition 


rf) 710 

P 


Parameter equal to 



PHIPTB 


SECTION II 


FORTRAN SYMBOLS 


FORTRAN 

SYMBOL NUMBER Definition 


AA 720 Gram atoms of element J in specie I. This 

two-dimensional array defines the chemical formula 
of the individual species. For example, if the 
order of the elements are assigned as follows: 

J = 1 Carbon 
J * 2 Hydrogen 
J = 3 Oxygen 
J = 4 Nitrogen 

Then, for this 4 element system (mm = 4) 
if the specie methane (CH^) is designated specie 
number 1. 

AA (1, 1) = 1.0 

AA (1, 2) = 4.0 

AA (1, 3) = 0. 

AA (1, 4) = 0. 

The resulting array for the species H, H , CH, , 
N 2 , NH 3 , NO, vI^O and CO is as follows: " "* 




C 

H 

0 

N 


CH. 

4 

1.0 

4.0 

0 . 

0 . 


H 

0 . 

1 . 

0 . 

0 . 


H 2 

0 . 

2. 

0 . 

0 . 


nh 3 

0 . 

3. 

0 . 

1 . 


NO 

0 . 

0 . 

1 . 

1 . 


h 2 o 

0 . 

2. 

1 . 

0 . 


CO 

I. 

0 . 

1 . 

0 . 

ABEXP 

722 

Pre-exponential factor for computing surface mass 
loss rate by sublimation 

ABUT 

724 

Initial value of 

the 

square root of pressure 
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FORTRAN 

SYMBOL 

NUMBER 

Definition 

AC0EF 

726 

Negative blocking coefficients when boundary layer 
flow through porous char is considered 

AEF 

728 

Activation energy of the forward reaction of 
pyrolysis gas kinetics 

AER 

730 

Activation energy of the reverse reaction of 
pyrolysis gas kinetics 

AF 

732 

Frequency factor of the foreward reaction of 
pyrolysis gas kinetics 

AI , BI , Cl : 
DI.KI.FI. 

, 733 

► 

Empirical constants for the enthalpy of the chemical 
specie of the pyrolysis gas above 1000 K by the 
equation AI*T+BI*T**2/ (2*1.987) + CI*T**3/(3*1.987) 

+ DI*t**4/ (4*1. 987) + EI*T**5/ (5*1.987) + 1.987*FI 

AIIjBII, 
CIIjDII, 
Eli, FII 

734 

Empirical constants for the enthalpy of the chemical 
specie of the pyrolysis gas below 1000 K as shown 
in AI, BI, Cl, DI, El, FI 

AR 

745 

Frequency factor of the reversed reaction of 
pyrolysis gas kinetics 

All — A55 

746 

Empirical constants of heat capacity of elements 
below 1000 K. 

BBEXP 

750 

Power of the exponential term for computing 
surface mass loss rate by sublimation 

CHANGE 

760 

A computing time saving device. Chemical reaction 
rates are not calculated until T(l) has 
changed more than CHANGE 

CMASSF 

770 

The mass fraction of the original material that 
is carbon 

C0N 

780 

Normalyzing constant for flow through porous media 

CS0LID 

790 

The mass fraction of a pyrolysis specie that is 
converted to char upon pyrolyzing 

DELH 

810 

Enthalpy of formation at 298.15 K 

DELWF0 

820 

Initial coolant flow rate 


70 



FORTRAN 

SYMBOL 

NUMBER 

Definition 

DUM1TB 

821 

Blank table (not used) 

O' 

w 

840 

Alphanumeric input of pyrolysis chemical reaction 
equation, for example, the reaction CH^ * C + 21^ 
would be input as CH4 = C + 2H2 

ERR0R 

850 

The test of T(l) for convergence from one itera- 
tion to the next. Convergence is assumed to have 
occurred when 



T W - ERR0R 

T(l) 



Where T(l) is from the previous iteration 

FRE0 

853 

Cyclic printing interval 

IDEBUG 

853 

When this input is greater than zero, extra 
printing will be done each iteration to aid debugging 

I0PT 

854 

Initial option for computing surface recession. 
There are 5 options available. These may be 
selected by assigning values to I0PT of 1, 2, 
3,4, or 5. When a value of 1 is assigned, 
surface recession is removed by (1) an 


exponential function of temperature or (2) by 
oxidation. Both forms are achieved by determining 
values for AEXP and BEXP . 

When a value of 2 is assigned, the mass loss 
rate is computed as a function of time 
from table WCD0TT.If NCD0T is a negative number, 
the table WCD0TT can be used to compute 
surface recession as a function of surface 
temperature . 

When a value of 3 is assigned, the thickness 
of the ablator is specified as a function of 
time in XTAB. When a value of 4 is assigned, 
surface removal is by sublimation and I0PT1=1. 
When a value of 5 is assigned, the surface 
recession is computed taking into account 
internal flow through porous media. (See 
ref. 9.) If option 5 (I0PT-5) is not used, 
the following inputs are not required: GAM1TB, 
GAM2TB, XRTAB, PHIPTB, PHIETB, MXR, NXR, 

TTBPHI, MPHI, NPHI, PERM, R, CON, and AC0EF. 


FORTRAN 

SYMBOL 

NUMBER 


Definition 

I0PT1 

855 

Option for computing sublimation I0PT1 » 1 for 
sublimation at the outer surface 

ICPBAR 

857 

Option for simplified treatment of pyrolysis gas 
kinetics. ICPBAR = 1. Enthalpy of pyrolysis 
gas must be supplied in CPBTAB 

IC0DE 

858 

Phase identification of pyrolysis gas components 
IC0DE = 0 - specie is a gas 
IC0DE =1 - specie is a solid 

IPRINT 

860 

Option for printing contents of NAMELIST entitled 
NAME1. When no value is input or when IPRINT * 0, 
NAME1 is printed. 

KTXDP 

870 

Option for variable density and conductivity of 
the third layer 

KTXDP =» 0 Density constant and equal to R0DP(1) 
KTXDP # 0. R0DP and CKNDPT are input as a function 
of position 

MALPHA 

880 

Order of interpolation 

for ALPHA table 

MALPH1 

890 

Order of interpolation 

for AL1 table 

MALPH2 

900 

Order of interpolation for AL2 table 

MAXIT 

910 

Maximum number of interations performed before time 
step reduced to one-half current value 

MCD0T 

920 

Order of interpolation 

in WCD0T table 

MCKNDP 

930 

Order of interpolation 

In KDP table 

MCKNP 

940 

Order of interpolation 

in KP table 

MCP 

950 

Order of interpolation 

for CP table 

MCPBAR 

960 

Order of interpolation 

for CPBAR table 

MGPDP 

970 

Order of interpolation 

for CPDP table 

MCP1 
MCP 2 

980 

981 

Order of interpolation 
Order of interpolation 

In CPPT1 table 
in CPPT2 table 

MCP 7 

987 

Order of interpolation 

in CPPT7 table 

MDUM1 

990 

Order of interpolation 

for DUM1 table 
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FORTRAN 

SYMBOL NUMBER 


Definition 



MHC 

MUE 

MHW 

MK 

MM 

MPHI 

MPRAT 

MQC 

MQR 

MQRAT 

MQS 

MS MIN 

MTAU 

MTB 

MX 

MXPRES 

MXR 

NALPHA 

NALPH1 

NALPH2 

MCCD0T 

NCD0T 

NCKNDP 


1000 Order of interpolation in HCTAB table 

1010 Order of interpolation in HETAB table 

1020 Order of interpolation in HW table 
1030 Order of interpolation in K table 
1040 Number of elements in pyrolysis gas species 
1050 Order of interpolation in PHIETB table 
1060 Order of interpolation in PRATT table 

1070 Order of interpolation in QC table 

1071 Order of interpolation in QR table 

1072 Order of interpolation in QRATT table 

1080 Order of interpolation in QS table 

1090 Minimum value of sublimation for which sublimation 
program STAT1 is used to decrease computing time 
when sublimation is occuring 

1100 Order of interpolation in TAUTAB table 

1110 Order of interpolation in TB table 

1120 Order of interpolation in X table 

1130 Order of interpolation in XPRESS table 

1140 Order of interpolation in XRTAB table 

1150 Number of entries in T ALP HA and ALPHA! 

1160 Number of entries in TTAL1 and AL1TAB 

1170 Number of entries in TTAL2 and AL2TAB 

1180 Number of entries in TTMCDT and WCD0TT. 

Negative value required by program logic 

1190 Number of entries in TTMBCT and NCD0TT 

1210 Number of entries in TCKNP and CKNDPT 
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FORTRAN 

SYMBOL 

NUMBER 

NCKNP 

1220 

NCP 

1230 

NCPBAR 

1240 

NCPDP 

1250 

NCP1 
NCP 2 

1260 

1261 

NCP7 

1267 

NDUM1 

1270 

NEQ 

1280 

NHC 

1300 

NHC0MB 

1310 

NPHC0MB 

1311 

NHE 

1320 

NHW 

1330 

NK 

1340 

NPEX 

1350 

NPHI 

1360 

NPYR 

1370 

NPRAT 

1380 

NQC 

1400 

NQR 

1410 

NQRAT 

1420 

NQS 

1430 

NREX 

1440 


Definition 

Number of entries in TTCKNP and CKNPTB 

Number of entries in TTABCP and CPTAB 

Number of entries in TTBCPB and CPBTAB 

Number of entries in TTCPDP and CPDPTB 

Number of entries in T1CPP and CPPT1 
Number of entries in T2CPP and CPPT2 

Number of entries in T7CPP and CPPT7 

Number of entries in vacant table 

Number of chemical reactions in pyrolysis gas kinetics 

Number of entries in TTABHC and HCTAB 

Number of entries in TTHC0MB and HC0MBTB 

Number of entries in PHC0MB 

Number of entries in TTABHE and HETAB 

Number of entries in TTABHW and HWTAB 

Number of entries in TTABK and CKTAB 

Power to which the concentration of the products 
is raised in pyrolysis gas kinetics reactions 

Number of entries in TTBPHI, PHIPTB and PHIETB 

Number of pyrolysis reactions 

Number of entries in TTPRAT and PRATT 

Number of entries in TTABQC and QCTAB 

Number of entries in TTABQR and QRTAB 

Number of entries in TTQRAT and QRATT 

Number of entries in TTABQS and QSTAB 

Power to which concentrations of reactants is 
raised in pyrolysis gas kinetics reactions 
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FORTRAN 

SYMBOL 

NUMBER 

Definition 

NS 

’390 

Number of solids in kinetic equations of pyrolysis 
gas 

NTAU 

1450 

Number of entries in TTAU and TAUTAB 

NTB 

1460 

Number of entries in TTABTB and TBTAB 

NTOPT 

1470 

New surface recession option changed to at T0PT 

NUMZ 

1480 

Number of Z values used 

NX 

1490 

Number of entries in TTABX and XTAB 

NX0PT 

1500 

New surface recession option changed to at X0PT 

NXPRES 

1510 

Number of entries in TTXPRE and XPREST 

NXR 

1520 

Number of entries in XRTAB, GAM1TB and GAM2TB 

0&DER 

1530 

ORDER 

Order of pyrolysis reaction m^ = (p) (AP) exp 

PELAM 

1550 

Pressure in ELAM table 

PERM 

1551 

Permeability of char 

PHC0MB 

1560 

Pressures in HC0MBTB 

PfiATT 

1570 

The ratio of the local pressure to the stagnation 
pressure 

PREHE 

1580 

Option for using punched input from program 
0PTUM (PREHE » 1.) 

PSC 

1581 

Stoichiometric coefficients of the products 
of pyrolysis gas kinetics equations. This array 
in conjunction with RSC (Stoichiometric coeffi- 


dents of the reactants) and SNAME (specie 
Identification and ordering) define the chemical 
reactants a~ they are used in calculation. The 
equations are written symbolically as: 


-BP/T) 


FORTRAN 

SYMBOL 


P5 

QRATT 

QSTAB 

R 

REFF 

RH0Z 

R0DP 

RSC 


NUMBER 


Definition 


S 


S 

N 

• rpcr 

N 

A 

r c J 

A 

M 


M 

E 


E 


to usi a simple example, if a system of only 4 
specie? are considered: H, H 2 , CH, and C 

The eqi ations 

CH. * C + 2 H. 

4 * 2 

H„ * 2H 

2 

are writtt a as: 


0 0 10' 

H 

*0 2 0 1' 

H 

,0100 

H 2 

" 2 0 0 0 

h 2 


CH. 

4 


CH. 

4 


C 

' 

C 


1585 Pressure tabl * for HCTAB 

1620 The ratio of 1 ocal convective heating rate to stag- 
nation convective heating rate. If QRATT set to zero, 
the pressure 1 computed using QS from QSTAB as 
follows: P w « R1 ?F*PC0EF*PRAT*QS/HE 

1625 Stagnation heati lg rate table used to obtain 
hypersonic press? re when QRATT-0. 

1630 Maximum radius of xblating body. Required only 
when flow through j. irons media effects are 
considered (I0PT ■ . 'I 

1700 Effective nose radius 

1640 Initial density of each component of material that 
will pyrolyze 

1641 Density of insulation layer 

1650 Stoichiometric coefficients of the reactants of 

pyrolysis gas kinetics equations See explanation of PSC 
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Definition 

SF 

1660 

Power on the temperature In the forward rate 
constant in pyrolysis gas reaction 

SNAME 

1670 

Alphanumeric input of specie names in pyrolysis 
gas kinetics. See explanation of PSC 

SR 

1680 

Power on the temperature in the reversed rate 
constants in pyrolysis gas reactions 

S1-- S5 

1690 

Empirical constants for heat capacity of 
elements below 1000 K 

TALPHA 

1710 

Temperature table corresponding to values in ALPHA! 

TAU0 

1720 

Starting time 

TAUTAB 

1730 

Table of time steps 

TBTA3 

1740 

Temperature to which back surface is radiating 

TCHAN 

1750 

Temperature change required before pyrolysis gas 
composition is recomputed 

TCKNDP 

1760 

Temperature table for CKNDPT 

TNBAR 

1770 

Back surface temperature at which coolant 
becomes active (TNST0P * 0). When TNST0P m 1, 
TNBAR is back surface temperature at which program 
will stop 

TNST0P 

1780 

Controls point at which computation stops. If TNST0P *0, 

program stops at ENDTAU. If TNST0P « 1, 

program stops when back surface temperature reaches TNBAR 

T0PT 

1790 

Time at which surface removal option may be 
changed using NT0PT 

TTABCP 

1800 

Temperature table for CPTAB 

TTABHE 

1810 

Time table for HETAB 

TTABHC 

1820 

Temperature table for HCTAB 

TTABHW 

1830 

Temperature table for HWTAB 

TTABK 

1840 

Temperature table for CKTAB 

TTABQC 

1850 

Time table for QCTAB 
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Definition 

TTABQR 

1860 

Time table for QRTAB 

TTABQS 

1870 

Time table for QSTAB 

TTABTB 

1880 

Time table for TBTAB 

TTABX 

1900 

Time table for XTAB 

TTAL1 

1910 

Time table for AilTAB 

TTAB2 

1920 

Time table for AL2TAB 

TTAU 

1930 

Time table, for TAUTAB 

T1BCPB 

1940 

Temperature table for CPPTAB 

T1CPP 

l*5u 

Temperature table for CPPT1 

T2CPP 

1951 

Temperature table for CPPT2 

T7CPP 

1957 

Temperature table for CPPT7 

TTBPH1 

1960 

Temperature table for PHIPTB and PHIETB 

TTCPDP 

1970 

Temperature table for CKNDPT 

TTCKNP 

1980 

Temperature table for CKNPTB 

TTDUM1 

1990 

Blank table 

TTELAM 

2000 

Temperature table for ELAMTB 

TTHC0MB 

2010 

Temperature table for HC0MBTB 

TTMCDT 

2020 

Time table for WCD0TT 

TTPRAT 

2030 

Time table for PRATT 

TTQRAT 

2040 

Time table for QRATT 

TTXPRE 

2050 

Time table for XPREST 

TZ 

2060 

Temperature at fixed locations (Z) 

WCCD00 

2070 

Initial value of WCCD0T 

WCD0T 0 

2080 

Initial value of WCD0T 
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WCD0TT 

2090 

Table used with option 2 

WCSD00 

2100 

Initial value of WCSD0T 

WCSD0T 

2101 

Surface mass loss rate by sublimation 

WPD0T0 

2110 

Initial value of WPD0T 

WPSUM 

2120 

Total mass flow past a given poini (one value for 
each station) 

XFW 

2141 

Molecular weight of each specie in pyrolysis gas 

XKIN 

2150 

X<XMIN, program stops 

X0PT 

2160 

Thickness of X at which surface removal option 
will change using NX0PT 

X0RDER 

2170 

Order of oxidation reaction 

XPMIN 

2180 

XP<XPMIN, program stops 

XPREST 

2190 

Square root of the total stream pressure in 
atmospheres, when considering subsonic or super- 
sonic flow. Set to .'.eso for comj mted pressure 

XRTAB 

2200 

Table to input radius of curvature as a func- 
tion of time when option 5 is used. 

XTAB 

2210 

Table to input char thickness as a function of 
time when option 3 is used 

YDEL 

2220 

Stefan Boltzmann constant in appro- 
priate units. 

YI 

2230 

The mole fraction of gas species produced by 
each pyrolyzable component 

Z 

2240 

Location of fixed stations measured from back 
surface 


TABLE II 


This table contains a list of input and output FORTRAN variables 
in the order they appear in the I/O lists. Instructions for use of library 
subroutines related to I/O are included in Appendix A. All tables use 
subroutine FTLUP for table look-up unless otherwise noted below under 
comments. 

The following input3 are entered by namelist entitled NAME1: 


FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 

NAME IN 

TABLE I 

SYMBOL 

INPUTS REQ'D 


(10) T 

460 

T 

I+J+M 

Used to initialize 
temperature distribution 

(20) Z 

580 

y-yo 

NUMZ 


(30) TZ 

2060 

T 

NUMZ 

Used to initialize tempera 
tures at fixed locations 
within solid (e.g. thermo- 
couple locations) 

(40) MCP 

950 

NA 

1 


(50) NCP 

1230 

NA 

1 

Table to input specific 

(60) TTABCP 

1800 

T 

NCP 

heat of char 

(70) CPTAB 

91 

C 

P 

NCP 


(80) MCP1 

980 

NA 

1 


(81) MCP2 

981 

NA 

1 


__ 


__ 

_ 


(86) MCP 7 

986 

NA 

1 


(90) NCP1 

1260 

NA 

1 

Table to input specific 

(91) NCP 2 

1261 

NA 

1 

heat of each componert in 
the virgin state 



— 

— 

- 


(96) NCP 7 

1266 

NA 

1 


(100) T1CPP 

1950 

T 

NCP1 


(101) T2CPP 

1951 

T 

NCP 2 



— — 

— 

— 


(106) T7CPP 

1956 

T 

NCP 7 
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FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENT'S 

NAME IN TABLE I 

SYMBOL 

INPUTS REQ'D 

(110) CPPT1 100 

C' 

P 

NCP1 


(111) CPPT2 101 

c , 

NCP2 


— —— 




(116) CPPT7 106 

V 

NCP7 


(120) MCPDP 970 

NA 

1 


(130) NCPDP 1250 

NA 

1 

Table to Input specific 

(140) TTCPDP 1970 

T 

NCPDP 

heat jf insulation layer 

(150) CPDPTB 110 

II 

C 

P 

NCPDP 


(160) MCPBAR 960 

NA 

1 


(170) NCPBAR 1240 

NA 

1 

Table to input enthalpy 
of pyrolysis gas as a 

(180) TTBCPB 1940 

T 

NCPBAR 

function of temp, when 
kinetic treatment of 

(190) CPBTAB 120 

C 

P 

NCPBAR 

pyrolysis gas is not 
calculated internally. 

This table not required 
when pyrolysis gas kinetics 
are treated in detail 
(ICPBAR^l) 

(200) MK 1030 

NA 

1 


(210) NK 1340 

NA 

1 

Table to input thermal con- 
ductivity of char as a func- 

(220) TTABK 1840 

T 

NK 

tion of tempera ture. 

(230) CKTAB 250 

k 

NK 


(240) MCKNP 940 

NA 

1 


(250) NCKNP 1220 

NA 

1 

Table to input composite 
thermal conductivity of 

(260) TTCKNP 1980 

T 

NCKNP 

virgin material as a func- 
tion of temperature 

(270) CKNPTB 260 

k' 

NCKNP 

(280 MCKNDP 930 

NA 

1 

Table input thermal 

conductivity of insulation 

(290) NCKNDP 1210 

NA 

1 

layer as a function of 
temperature and position 
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FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 


NAME 

IN 

TABLE I 

SYMBOL 

INPUTS REQ’D 



(300) 

TCKNDP 

1760 

T 

NCKNDP 



(310) 

CKNOPT 

270 

k" 

NCKNDP 



(320) 

MHW 

1020 

NA 

1 

Table to input enthalpy of 
external gas stream at the 


(330) 

NHW 

1330 

NA 

1 

surface temperature as a 
function of temperature 


(340) 

TTABHW 

1830 

T 

N>iW 



(350) 

HWTAB 

210 

H 

w 

NHW 



(360) 

MALPHA 

880 

NA 

1 

Table to input absorptivity 


(370) 

NALPHA 

1150 

NA 

1 

of the surface as a 
function of temperature 


(380) 

TALPHA 

1710 

T 

NALPHA 



U90) 

ALPHAT 

590 

“t 

NALPHA 



(400) 

MALPH1 

890 

NA 

1 

Tabic to input blocking 
factor of surface reaction 


(410) 

NALPH1 

1160 

NA 

1 

products as a function of 
time 


(420) 

TTAL1 

1910 

T 

NALPH1 



(430) 

AL1TAB 

(i00 

c 

NALPH1 



(440) 

MALPH2 

900 

NA 

1 

Table to input blocking 
factor of pyrolysis reaction 


(450) 

NALPH2 

1170 

NA 

1 

products as a function of 
time 


(460) 

TTAL2 

1920 

T 

NALPH2 



(470) 

AL2TAB 

610 

a 

P 

NALPH2 



(480) 

MHE 

1010 

NA 

1 

Table to input enthalpy of 

i 

NHE 

1320 

NA 

l 

gas external to boundary 
layer as a function of 
time 

(490) 


(500) 

TTABHE 

1810 

t 

NHE 


H 

(510) 

HETAB 

200 

H 

e 

NHE 


: 

(520) 

MQC 

1070 

NA 

1 

Table to input convective 
heating rate as a function 

j 

(530) 

NQC 

1400 

NA 

1 

of time 


32 




i 




1 

i 

5 


4 


J* 

* 


i 

i 



FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 

NAME 

IN 

TABLE I 

SYMBOL 

INPUTS REQ'D 


(540) 

TTABQC 

1850 

t 

NQC 


(550) 

QCTAB 

390 

q c 

NQC 


(560) 

MQR 

1071 

NA 

1 

Table to input radiant 
heating rate to the surface 

(570) 

NQR 

1410 

NA 

1 

as a function of time 

(580) 

TTABQR 

1360 

t 

NQR 


(590) 

QRTAB 

400 

q r 

NQR 


(600) 

MCD0T 

920 


1 

Table of surface recess- 
ion rate as a function of 

(610) 

NCD0T 

1190 


1 

time when option 2 is used 
6 NCCD0T is greater than 

(620) 

TTMCDT 

2020 

T 

NCD0T 

or equal to zero. When 




O 


NCCD0T is less than zero. 

(630) 

WCD0TT 

2090 

m 

NCD0T 

this table is used to 






input surface recession as 
a function of surface 
temperature. 

(640) 

MX 

1120 

NA 

1 

Table to input char thick- 
ness as a function of time 

(650) 

NX 

1490 

NA 

1 - 

when option 3 is used 

(660) 

TTABX 

1900 

t 

NX 


(670) 

XTAB 

2210 

X 

NX 


(680) 

MDUM1 

990 


1 


(690) 

NDUM1 

127 


1 

Table for expansion 
(not used) 

(700) 

TTDUM1 

1990 


NDUM1 


(710) 

DUM1TE 

821 


NDUM1 


(720) 

MQS 

1080 

NA 

1 

Table of stagnation heating 
rate as a function time 

(730) 

NQS 

1430 

NA 

1 


(740) 

TTABQS 

1870 

T 

NQS 


(750) 

QSTAB 

410 

q s 

NQS 
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FORTRAN ITEM NUMBER 
NAME IN TABLE I 

PHYSICAL 

SYMBOL 

NO. OF 
INPUTS REQ'D 

COMMENTS 

(760) FI 

160 

.F 

1 



(770) R0 

— 




Calculated Internally 

(780) ROP 

— 




Calculated internally 

(790) RQDP 

680 

P" 

1 


can have 10 vilues which are 
specified for 10 equally 
spaced thicknesses of 
insulation numbered from the 
back surface of the ablator. 
Either one value or 10 values 
required. 

(800) DELHF 

240 

AH f 

1 



(810) DELHP 

— - 

AH 

P 

1 


Calculated internally 

(820) EIPJ 

640 

E i+j 

1 



(830) CIPJ 

130 

C i+j 

1 



(840) EIJM 

650 


1 



(850) CIJM 

140 

C i+j+m 

1 



(860) BETA 

620 

8 

1 



(870) ETA 

660 

n 

1 



(880 % EPSONE 

630 

e i 

1 



(890) TNST0P 

1780 

NA 

1 



(900) X 

540 

X 

1 



(910) MHC 

1000 

NA 

0 


Not used for input. 

(920) NHC 

1300 

NA 

0 



(930) TTABHC 

1820 

T 

5 


A value of 5 is assigned 
internally. Not used for input. 

(935) P5 

1585 


5 


A value of 5 is assigned 


Internally. Mot used for Input. 
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COMMENTS 


FORTRAN ITEM NUMBER PHYSICAL NO. OF 

NAME IN TABLE I SYMBOL INPUTS REQ'D 


(940) HCTAB 190 He 25 


(950) MIB 

1110 

FA 

1 

(960) NTB 

1460 

NA 

1 

(970) TTABTB 

1880 

t 

NTB 

(980) TBTAB 

470 

t b 

NTB 

(990) XP 

550 

X' 

1 

(1000) XDP 

560 

X" 

1 

(1010) TNBAR 

17 VO 

T 

i+j 

or 

T 

i+j+rn 

1 

(1020) I 

241 

i 

1 

(1030) J 

242 

j 

1 

(1040) M 

300 

m 

1 

(1050) TAU0 

490 


1 

(1060) WCD0T0 2080 


1 

(1070) WPD0T0 2110 


1 

(1080) DELWF0 

820 


1 

(1090) FRE0 

851 


1 

(1100) ENDTAU 

500 

fc f 

1 

(1110) X0PT 

2160 


1 

(1120) I0PT 

854 


1 


Not used for input. A 
value of 25 is assigned 
internally. Table to input 
latent heat of surface removal 
process as a function of 
temperature and 
pressure. Uses subroutine 
DISCOT for table look-up. 


Table to input tempera- 
ture to which back surface 
is radiating as a func- 
tion of time 
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FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 


NAME IN TABLE I 

SYMBOL 

INPUTS REQ'D 


- 

(1130) NX0PT 

1500 


1 



(1140) XMIN 

2150 


1 



(1150) XPMIN 

2180 


1 



(1160) NUMZ 

1480 


1 



(1170) T0PT 

1790 


1 



(1180) NT0PT 

1470 


■1 



(1190) AEXP 

10 


1 



(1200) BEXP 

40 


1 



(1210) APEXP 



0 

Do not use 


(1220) BPEXP 



0 

Do not use 


(1230) X0RDER 

330 


1 



(1240) CE 

90 

C 

e 

1 


• 

(1250) ELAM 

— 


o 

Do not use 


(1260) HC0MB 

— 


l 

Do not use 


(1270) REFF 

1700 





(1280) PC0EF 

180 

VT 

l 

square root of constant 
in hypersonic pressure term 
A value of 7.5 in SI units 
and 20.27 in English units 
is recomnended. 




(1290) MPRAT 

1060 

NA 

l 

Table to input ratio of 
local pressure to stag- 


(1300) NPRAT 

1380 

NA 

l 

nation pressure as a 
function of time 


(1310) TTPRAT 

2030 

t 

NPRAT 



(1320) PRATT 

1570 

P/P 

s 

NPRAT 
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FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 

NAME IN TABLE I 

SYMBOL 

INPUTS REQ'D 


(1330) MXPRES 1130 

NA 

1 

Table to input square 
root of the total stream 

(1340) NXPRES 1510 

NA 

1 

pressure as a function 
of time 

(1350) TTXPRE 2050 

t 

NXPRES 


(1360) XPREST 2190 

Vp 7 

s 

NXPRES 


(1370) MQRAT 1072 

NA 

1 

Table to input the 
ratio of local laminar 

(1380) NQRAT 1420 

NA 

1 

to stagnation convec- 
tive heating as a func- 

(1390) TTQRAT 2040 

■ t 

NQRAT 

tion of time 

(1400) QRATT 1620 

q/q s 

NQRAT 


(1410) WCSD00 2100 


1 


(1420) WCCD00 2070 


1 


(1430) AL 280 

1 

1 


(1440) MTAU 1100 


0 

Make no entry 

(1450) NTAU 1450 

NA 

1 

Table co input time 
step as a function of 

(1460) TTAU 1930 

t 

NTAU 

time with zero order 
interpolation (time step 

(1470) TAUTAB 1730 

At 

NTAU 

remains constant at TAUTAB 
value until time equals or 
exceeds next TTAU value.) 
All entries in this table 
should be powers of 2, e.g 
1/32, 1/16, 1/2, 1, 2, etc 

(1480) TCHAN 1750 


1 

A value of 25 is 
recommended . 

(1490) MAXIT 910 


1 

A value of 3 is 
recommended 

(1500) ERR0F 850 


1 


(1510) IDEBUG 853 


1 


(1520) THETA — 


0 

Co not use 

(1530) ABEXP 722 


1 
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FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 

name in TABLE I 

SYMBOL 

INPUTS REQ'D 


(1535) BBEXP 750 


1 


(1540) I0PT1 855 


1 


(1550) PREHE 1580 


1 


(1560) ABUT 724 


1 


(1570) MXR 1140 

NA 

1 

Table to input approxi- 
mate geometric coefficients 

(1580) NXR 1520 

NA 

1 

as a function of position 
(see reference 9) 

(1590) XRTAB 2200 

x/r 

NXR 


(1600) GAM1TB 621 

X 

NXR 

Uses (1570 ) , (1580) , and 
(1590) above to input 

(1610) GAM2TB 621 

^2 

NXR 

approximate geometric 
coefficient as a function 

(1620) NPHI 1360 



of position 

Table to input parameter 
as a function of tempera- 

(1630) MPHI 1050 



ture (see reference 9) 

(1640) TTBPHI I960 


NPHI 


(1650) PHIPTB 710 

v 

NPHI 


(1660) PHIETB 700 

M 

V 

NPHI 

Use (1620), (1630), and 
(1640) above to input 
pnrameter 99 a function 
of temperature 



(1670) PERK 1551 

k 

1 

See reference 9 

(1680) R 420 



See reference 9 

(1690) C0N 780 



See reference 9 

(1700) AC0EF 726 


3 


(1710> YDEL 2220 


1 

. 0.6 



Table to input aN„ as a 

(1720) ELAKTB 350 

X 

40 

function of surface 
temp, and pressure 

(1730) ITFlAh 2000 

T 

10 

Uses D1SC0T for 
table look-up 
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FORTRAN ITEM NUMBER 
NAME IN TABLE I 

PHYSICAL 

SYMBOL 

NO. OF 
INPUTS REQ'D 

(1740) PELAM 1550 

P 

4 

(1750) IPRINT 860 



(1760) KTXDP 870 



(1770) ICPBAR 857 



(1780) MSMIN 1090 



(179O)TTHC0MB 2010 

T 

28 

(18OO)HC0MBTB 841 

AH . 
comb 

7 

(1810) PHC0MB 1560 

P 

4 

(1820) NHC0MB 1310 

NA 

0 

(183O)iiPHC0MB 1311 

NA 

0 

(1840) CMASSF 770 


1 

(1850) CHANGE 760 


1 

(1860) DDELHP 220 

AH 

P 

NPYR 

(1870) CS0LID 790 


NPYR 

(1880) RH0Z 1640 

P z 

NPYR 

(1890) NPYR 1370 


1 

(1900) AP 20 


NPYR 

(1910) BP 50 


NPYR 

(1920) RH0 665 

P 


(1930) RH0C 670 

p c 

I+J 

(1940) 0RDER 1530 


NPYR 

(1950) WFSUM 2120 

m 

P 

I+J 


COMMENTS 


Table to input heat of combus- 
tion of the char 
as a function 
of temperature and 
pressure. Uses DISCOT 
for table look-up. 

Set internally to a 
value of 28 

Set internally to a 
value of 4 
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FORTRAN 

NAME 


ITEM NUMBER 
IN TABLE I 


PHYSICAL 

SYMBOL 


NO. OF 
INPUTS REQ'D 


COMMENTS 


THE FOLLOWING INPUTS ARE ENTERED BY INDIVIDUAL FORMAT IN SUBROUTINE IN0UT 


(1960) 

NG 

340 

1 X 


(1970) 

NS 

1390 

1 

> F0RMAT (416) 

(1980) 

MM 

1040 

1 

(1990) 

NEQ 

1280 

1 


(2000) 

EQN 

840 

12 

F0RMAT (12A6) 

(2010) 

RSC 

1650 

Variable 

FORMAT (6F4.0) input 
required for each term 

(2020) 

PSC 

1581 

Variable 

in chemical reaction 


e.g. for the cheuical 

reaction CH. * C, v + 
4 (s) 

2H 2 , 1 RSC and 2 PSC 

are required 


(2030) 

NREX 

1440 

Same 


F0RMAT (514.0) 




as PSC 


- 

(2040) 

NPEX 

1350 

Same as 






PSC 



(2050) 

AF 

732 

1 



(2060) 

SF 

1660 

1 


> P0RMAT (E8.0,1X,2F6.0) 

(2070) 

AEF 

728 

1 

j 


(2080) 

AR 

745 

1 

' 


(2090) 

SR 

1680 

1 


) F0RMAT (E8.0.1X, 2F6.0) 

(2100) 

AER 

730 

1 



REPEAT 

THE SEQUENCE OF ITEMS 

2000 THROUGH 2100 

* 

FOR EACH CHEMICAL REACTION 

(2110; 

SI 

1690 

1 

\ 







FORMAT (7E10.5) One set 

(2120) 

S2 

1690 

1 


of inputs for each 






element in the system 

(2130) 

S3 

1690 

1 


) (l.e. MM cards) 

(2140) 

S4 

1690 

1 
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FORTRAN ITEM NUMBER 

PHYSICAL 

NO. OF 

COMMENTS 

NAME IN 

TABLE I 

SYMBOL 

INPUTS REQ’D 


(2150) S5 

1690 


X J 


(2160) All 

746 


1 \ 

FORMAT (7E10.5) One set 

(2170) A22 

746 


1 

of inputs for each 
element in the system 

(2180) A33 

746 


1 

> (i.e. Ml cards) 

(2190) A44 

746 


1 


(2200) A55 

746 


1 J 


(2210) XFW 

290 


1 

FORMAT (20X, E10.4, 10X 

(2220) SNAME 

1670 


1 

A6, 16, E15.4) 

^ KSPI sets of input 

(2230) IC0DE 

858 


1 

(2240) DELH 

810 


1 ' 


(2250) YI 

2230 


NPYR*KSP1 

FORMAT (7E10.4) 

KSPI cards with NPYR 
inputs per card 

(2260) AA 

720 


MM*KSP1 

FORMAT (7E10.5) 

(2270) AI 

733 


1 ] 


(2280) BI 

733 


i 


(2290) Cl 

733 


1 


(2300) DI 

733 


1 


(2310) El 

733 


1 

) FORMAT (7E10.5) 

(2320) FI 

733 


1 


(2330) GI 

733 


1 


(2340) All 

734 


1 


(2350) BII 

734 


1 


(2360) CII 

734 


i 

> FORMAT (7E10.5) 

(2370) DII 

734 


1 
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FORTRAN ITEM NUMBER 
NAME IN TABLE I 


PHYSICAL NO. OF 
SYMBOL INPUTS REQ'D 


COMMENTS 


(2380) Eli 734 
(2390) FII 734 
(2400) Gil 734 


1 
1 

1 / 


THE SEQUENCE OF ITEMS 2260 THROUGH 2400 IS REPEATED KSP1 TIMES 
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Preliminary Printing Before Computation Begins 

Output from NAMELIST entitled NAME1: This list is printed 
automatically unless IPRINT is different from zero. The list is headed 
by $NAME1 and contains each entry in NAME1 followed by the actual data 
contained in the locations reserved f^r that variab]e. If all locations 
reserved are not input, the remaining ones will be filled with zeros. 
This list is followed by $END. 

Title card: The contents of the title card are printed next and 

followed by the title "INPUT DATA." Some selected tables with titles 
are printed. These tables are followed by some items that remain 
constant during computation. 

Cyclic Printing 

Printing interval ; Cyclic printing of values that change with 
time are printed at intervals of FRE0. 

Items printed each cycle ; The number of iterations needed to 
converge during the last time step is printed last. The current time 
TAU and the last time step DELTAU used in the computation are printed 
first. These items are followed by 10 columns of numbers. The first 
column on the left is the temperature of each finite difference station 
from the surface inward. The second column (WWPD0T) is the mass flow 
rate of pyrolysis gas past each station. The third column RH0C is the 
char density at each station. The next seven columns are the local 
density of each of the solid components that make up the original 
ablation material (some columns may have all zero entries for a simple 
material). 

These 10 columns are followed by a list of current values of 
variables that are individually titled. 

Additional printing (1) when IDEBUG i 0., there is additional 
P r ^ n *-^ n 8 that is done to aid debugging. The following items are printed 
each iteration: TAU, WCD0T , WCCD0T, WCSD0T, WPD0T, X, XP, QCNET, 
1EMPERATURES of each station during current and previous iterations. 

(2) When ICPBAR / 1, and the pyrolysis reactions are treated in detail, 
the following items are nrinted each interval of FRE0. A group of items 
dealing with front surface conditions. This group is followed by a list 
en fifl e d ENTHALPY, This list contains the enthalpy of the pyrolysis gas 
at each c inite inference station where gas flow is large enough to be 
significant. The next list contains information about each reacting 
specie in the pyrolysis gas. The final list in this group is entitled 
WDEP and gives the current rate of deposition of carbon from the gas 
phase that occurs at each finite-difference station. 


TABLE III 


alphabetical listing of fortran symbols 


FORTRAN 

SYMBOL 

PHYSICAL 

SYMBOL 

LOCATION 
IN TABLE I 

LOCATION 
IN TABLE II 

AA 


720 

2260 

AAL 




ABEXP 


722 

1530 

ABIT 




ABUT 


724 

1560 

AC0EF 


726 

1700 

AC0N1 




AC0N4 




AC0N9 




AEF 


728 

2070 

AER 


730 

2100 

AEXP 

A 

10 

1190 

AF 


732 

2050 

AI 


733 

2270 

All 


734 

2340 

AL 

1 

280 

1430 

ALAM 

> 



ALPHA 

a 

590 


ALPHAS 




ALPHAT 


590 

390 

ALPHA! 





ALPHA?. 



I 

FORTRAN 

PHYSICAL 

LOCATION 

LOCATION 

SYMBOL 

SYMBOL 

IN TABLE I 

IN TABLE 11 

AL1TAB 

c 

600 

430 

AL2TAB 

a 

P 

610 

470 

AMUT 

AMI 




AP 

A’ 

20 

1900 

AR 


745 

2080 

ASQ 

APEXP 



12x : 

AVGX* 

AZ 

AZETA 




A1 

All 


745 

2160 

A22 


746 

2170 

A33 


746 

2180 

A44 


746 

2190 

A55 

A2 


746 

2200 

B 

BAT 

BB 

BC 




BBEXP 


750 

1535 


95 



FORTRAN 

SYMBOL 

PHYSICAL 

SYMBOL 

LOCATION 
IN TABLE I 

LOCATION 
IN TABLE II 

BETA 

3 

620 

860 

BEXP 

B 

40 

1200 

BIT 


734 

2350 


BL0CK 


BL0CK1 

BN 

B0X 


BP 

B' 

50 

1910 

BP EXP 



1220 

C 


80 


CAT 




CC 




CE 

Ce 

90 

1240 

CELL 




Cr-LL1 




CHANGE 


760 

1850 

Cl 


733 

2290 

CII 


734 

2360 

CIJM 

C i+j+m 

140 

850 

CIPJ 

C. , . 
i+J 

130 

830 

CK 




CKDP 




CKNDPT 

k" 

270 

310 

CKNPTB 

k' 

260 

270 
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C-> 


FORTRAN 

PHYSICAL 

LOCATION 

LOCATION 

SYMBOL 

SYMBOL 

IN TABLE I 

IN TABLE II 

CKP 

CKTAB 

k 

250 

230 

CMASSF 


770 

1840 

C0LL 

C0N 

C0NST 

CP 

CP BAR 
CPBTAB 

C 

P 

120 

190 

CPB1 

CPB2 

CPDP 

CPDT 

CPDT1 

CPDPTB 

V' 

110 

150 

CPP 

CPPT1 

t 

c 

100 

110 

CPPT2 
CPPT7 
CP TAB 

p 

c 

p 

91 

70 

CS0LXD 


790 

1870 


CUE 

CZ 

D 
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FORTRAN PHYSICAL LOCATION LOCATION 

SYMBOL SYMBOL IN TABLE I IN TABLE II 

DABC 

DABUT 

DB2KSQ 

DBL0CK 

DCALL 

DCELL1 

DCOLL 

DC 

DDD 


DDELHP 


220 

1860 

DELH 


810 

2240 

DELHF 


240 

800 

DELHP 



810 

DELSUM 



• 

DELTAU 

At 

690 


DELWF 

Aw f 

530 


DELWFO 


820 
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Figure 2.- Schematic of finite difference network. 





